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ABSTRACT 


The  open-end  resonator  in  suspended  substrate  transmission  line  was 
analyzed  in  terms  of  two  open-end  discontinuities  on  a  line  segment.  A 
procedure  for  computing  the  resonant  frequency  and  constant  of  propa¬ 
gation  was  demonstrated  using  a  full  wave  analysis.  The  fringing 
capacitance  was  computed  using  an  equivalent  length  extension  model 
and  a  transmission  line  circuit  model.  The  characteristic  equation  was 
derived  using  Galerkin's  method  applied  in  the  Fourier  transform  do¬ 
main.  The  calculation  has  been  carried  out  in  two  increasing  orders  of 
approximation,  and  the  results  compared.  The  discontinuity  capacitance 
at  the  resonator  open  ends  was  calculated  for  a  range  of  line  dimensions 
and  substrate  dielectric  constants.  The  dispersive  propagation  constant 
of  suspended  substrate  line  was  also  calculated  by  the  Galerkin  method. 
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I.  INTRODUCTION 


The  analysis  of  planar  transmission  lines  in  the  Fourier  transform 
domain  or  spectral  domain  is  superior  to  many  numerical  methods  in  the 
spatial  domain.  The  analysis  in  the  Fourier  transform  domain  was  used 
earlier  by  Yamashita  and  Mittra  [Ref.  4]  for  computation  of  the  charac¬ 
teristic  impedance  and  phase  velocity  of  a  microstrip  line  based  on  a 
quasi-TEM  approximation.  A  variational  method  was  used  in  the  Fourier 
transform  domain  to  calculate  the  line  capacitance  from  the  assumed 
charge  density.  This  is  a  low-frequency  approximation  neglecting  longi¬ 
tudinal  electric  and  magnetic  fields  supported  by  the  strip  transmission 
line  having  two  dielectric  media  in  its  cross  section,  or  inhomogeneous 
line. 

As  the  operating  frequency  is  increased,  dispersion  characteristics  of 
the  inhomogeneous  line  become  important  for  precise  designs.  This  re¬ 
quirement  has  led  to  the  full  wave  analysis  of  microstrip  lines,  represented 
by  the  work  of  Denlinger  [Ref.  15],  who  solved  the  integral  equations 
using  a  Fourier  transform  technique.  The  solution  by  his  method,  how¬ 
ever,  strongly  depends  on  the  assrmed  current  distributions  on  the  strip 
in  process  of  solution.  To  avoid  this  difficulty  and  permit  systematic  im¬ 
provement  of  the  solution  for  the  current  components  to  a  desired  degree 
of  accuracy,  a  new  method  was  presented  by  Itoh  and  Mittra  and  com¬ 
monly  called  the  Spectral  Domain  Approach  (  SDA  ).  In  SDA, 
Galerkin's  method  is  used  to  yield  a  homogeneous  system  of  equations 
to  determine  the  propagation  constant  and  the  characteristic  frequencies 
from  which  the  equivalent  circuit  is  derived. 
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In  each  of  these  methods  the  Fourier  transform  is  taken  along  the 
coordinate  axes  in  the  plane  of  the  strip.  By  virtue  of  the  Fourier  trans¬ 
form  domain  analysis  and  Galerkin's  method,  SDA  has  several  features: 

•  Easy  formulation  in  the  form  of  a  pair  of  algebraic  equations. 

•  Variational  nature  in  determination  of  the  propagation  constant. 

The  Spectral  Domain  Approach  is  applicable  to  the  following  struc¬ 
tures: 

•  Most  planar  transmission  lines  such  as  mic'ostrip,  finline,  and 
coplanar  waveguide  in  multilayer  configurations. 

•  Both  open  and  enclosed  structures. 

•  Slow-wave  lines  with  lossy  dielectric  materials. 

•  Resonators  of  planar  configurations. 

An  efficient  formulation  of  the  SDA  is  achieved  in  the  present  work 
by  enclosing  the  suspended  substrate  in  a  metallic  shield  enclosure.  In 
this  way  the  Fourier  transform  becomes  the  finite  Fourier  transform, 
calculated  as  a  summation.  Singularities  of  the  integrand  can  be  avoided, 
and  the  truncated  summations  have  acceptable  convergence  behavior. 
As  applied  here,  the  analysis  assumes  loss-free  conductors  and  dielectric 
media,  and  also  assumes  infinitesimal  thickness  for  the  strip  conductors. 
The  resonant  frequencies  are  obtained  by  numerically  solving  the  char¬ 
acteristic  equation,  and  equivalent  open-circuit  capacitances  are  obtained 
using  the  end-effect  and  the  transmission-line  models.  The  details  of  the 
analysis  method  will  appear  in  Chapter  II  of  this  thesis. 

In  Chapter  III,  discussions  are  presented  concerning  the  accuracy  of 
the  solution  for  the  resonant  frequencies,  the  propagation  constants,  and 
the  end-effect  at  th?  open  end  of  the  microstrip  structure  by  calculation 
of  the  fringing  capacitance.  These  solutions  for  various  orders  of  ap- 
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proximation  are  compared,  based  on  the  assumptions  for  various  current 
distributions  in  x  and  z  -  directions. 
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II.  SPECTRAL  DOMAIN  APPROACH  FOR  SHIELDED 
SUSPENDED  SUBSTRATE  RESONATOR 

A.  INTRODUCTION 

The  shielded  suspended  substrate  resonator  to  be  analyzed  is  shown 
in  Figure  1.  A  rectangular  strip  of  width  w  and  length  1  is  placed  on  the 
suspended  substrate  which  is  separated  from  the  ground  plane  by  an  air 
gap.  The  sides  and  the  top  of  the  structure  are  enclosed  with  metallic 
shielding  walls.  Thus  the  entire  structure  is  considered  to  be  the  sus¬ 
pended  microstrip  resonator  located  in  enclosed  partially  filled 
waveguide.  It  is  assumed  that  the  thickness  of  the  strip  is  negligible  and 
that  all  the  media  and  conductors  are  lossless.  The  shielding  waveguide 
has  dimensions  3.2  mm  X  1.575  mm  X  100.0  mm,  which  corresponds  to 
a  system  in  actual  practice.  For  simplicity,  the  strip  is  assumed  to  be 
symmetrically  located. 

The  operating  frequency  is  chosen  to  be  below  cutoff  of  the  shielding 
waveguide  partially  filled  with  substrate  material,  to  avoid  interaction 
with  waveguide-mode  resonances.  From  the  formula  in  [Ref.  10]  the 
cutoff  frequency  for  the  structure  in  Figure  1  is  shown  in  Table  1. 

Table  1.  THE  CUTOFF  FREQUENCY  OF  SUSPENDED  SUBSTRATE  LINE 


Dielectric  constant  (e,) 

Frequency  cutoff  (  GHZ  ) 

2.2  (  duroid  ) 

44.7574 

9.8  (  alumina  =  AliO^  ) 

43.3333 

Figure  1,  Shielded  suspended  substrate  resonator 


B.  METHOD  OF  SOLUTION  FOR  RESONANT  FREQUENCY 
1.  Field  Equations 

The  analysis  of  field  equations  on  strip  transmission  lines  of  the 
type  considered  here  has  been  carried  out  by  Uwano  and  Itoh  [Ref.  2]. 
These  authors  treated  the  problem  by  applying  Galerkin's  method  of 
moments  in  the  Fourier  transform  domain.  Their  work  is  summarized 
in  a  set  of  Green's  function  equations  as  follows 


^xr^z  ^xx^x  ~  ^x 
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Where  the  tilde  over  the  symbols  indicates  the  Fourier  transformed 
quantities.  Here,  and  are  the  transformed  z  and  x  components  of 
the  currents  on  the  strip  conductor.  and  are  the  components  of 
electric  field  tangential  to  the  substrate  surface.  ,  Z^  are  the 

Green's  impedance  functions,  and  are  defined  as  the  following  : 


Z,,  =  '  ,  irz,  +  cc'zj 

a^  +  r 

(Ic) 

z,,  =  -  [z,  -  zj 

a +r 

(1^ 

z  =  z 

^xz  ^zx 

(1^) 

'  2  +  phj 

a"  +  j5 

(ly) 

^  Cr2Cr3  +  CtiCty/y2lyy\  +  +  yj3/>’v-2 

(ig) 

~  >'72^^2  +  >’z3^^3 

•^/i  2 

yzf/z2^^l^^2  +  yzl>'z3^^1^^3  +  )'z2)'z3^^2C^3  +  >’r2 

(i«) 

Cq  =  cor/z  ^i/z 

(1/) 

O2  =  coth  fit 

(ly) 

Ct^  =  coth  yyi 

(l/c) 

2  2  ^  «2  ,2 
y/  =  a  +  ^  -  /c/ 

(1/) 
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/  2  2 
ki  ^O)  nfii 


_ 

yyl  yi 


_ 

^2/  Zi 


yi^jcozi 

zi=jai^l 

El 

f^l  =  /^rA^o 


(1/w) 


(1/z) 


(lo) 

(1/j) 

(W) 

(Ir) 

(1^) 


Where  subscripts  /  =  1,2,3  refer  to  the  corresponding  regions  1,2,3,  co  is 
the  operating  frequency,  and  £„  and  are  the  free-space  permittivity  and 
permeability,  respectively.  In  this  thesis,  the  Fourier  transform  is  carried 
out  in  the  bounded  region  interior  to  the  enclosed  resonator.  In  order 
to  examine  the  layered  structures  in  Figure  1  with  added  bounding  verti¬ 
cal  electric  and  magnetic  walls  at  a:  =  ±  1.6mm  ,  z  =  ±  50mm  ,  the  finite 
Fourier  transform  should  be  used  [Ref.  20].  There  we  have  the  two- 
dimensional  finite  Fourier  transform  pairs 


fl-6  f50  .. 

W=  Ax.y,z)e  '”^dxdz 


(2u) 
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{2b) 


My.  ^)  = 

“'>  0n, 

Where  a„  and  (i„  are  the  discrete  transform  variables  for  the  dominant 
mode  and  are  defined  by,  a„  =  {n  —  l/2)Tr/^  for  even,  —H^  odd(in  x) 
modes,  and  ^„  =  {m—  \l2)njb  for  even,  odd(in  z)  modes.  Note 
that  b  is  the  half  length  of  the  resonator  and  is  defined  to  be  5  times  the 
length  of  strip,  and  a  is  the  half  width  of  the  resonator.  From  (2b)  we 
have 


z) = yy  ?(«„,>■, 
1 1: 


(2c) 


Equations  (la,b)  contain  four  unknowns  J,  ,  ,  E^  ,  and  E^  .  Two 

unknowns  and  E^  ,  however,  can  be  eliminated  by  applying  Galerkin's 
method  in  the  spectral  domain.  The  first  step  is  to  expand  the  unknowns 
and  Jx  in  terms  of  known  basis  functions  and 


A'2 

iVl 

Jx{^nPn^  ~  ^  ^(Pxk{^nPn^ 
k=\ 

Where  Cj  and  are  unknown  coefficients.  The  basis  functions  must  be 
chosen  to  approximate  the  true  but  unknown  distributions  on  the  strip. 
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Si.ice  the  current  is  nonzero  only  on  the  strip,  therefore  each  basis  func¬ 
tion  must  be  chosen  so  that  it  is  nonzero  only  on  the  strip.  After  substi¬ 
tuting  (3)  into  (la,b)  and  taking  inner  products  of  the  resulting  equations 
with  the  basis  functions  J^i  and  for  different  values  of  i,  this  process 
yields  the  equations, 


iVl 

A'2 

z 

‘J  J 

^4‘’^4=o 

,/  =  1,2,3,..,A^2. 

(4u) 

k=l 

Al 

A'2 

z 

Ay  Cj  -b 

,/•=  1,2,3,.., A^l. 

(4^) 

^=1 

k=l 

Where  from  definition  of  inner  products  associated  with  the  Fourier 
transform  defined  by  (2),  equations  (4)  lead  to  a  set  of  homogeneous 
equations  forming  the  solution  matrix  for  the  Cj  and  •  The  matrix  el¬ 
ements  are 


-Koo  +«= 

+  00 
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+00 

E  E  (Sc) 

>^=—oorn——oo 

+00 

=  ‘/*  E  E  JxKPmtixxKfimVA'^nfim)  (Sd) 

The  right  hand  side  of  (la,b)  can  be  eliminated  in  Galerkin's  process  via 
the  application  of  Parseval's  relation,  because  the  current  J^ix)  ,  J^ix)  and 
the  field  components  E^{x,D  +  T)  ,  EJ^x,D  +  T)  are  nonzero  in  the  com¬ 
plementary  regions  of  x.  Note  that  the  solution  matrix  is  symmetrical. 
In  order  to  have  non-trivial  values  for  ,  dj  ,  the  determinant  of  the 
matrix  must  be  zero.  From  this  we  have  the  solution  for  resonant  fre¬ 
quency.  The  accuracy  of  the  solution  can  be  systematically  improved  by 
increasing  the  number  of  basis  function  (N 1  +  N2)  and  by  solving  a  larger 
size  matrix  equation.  However,  if  the  first  few  basis  functions  are  chosen 
so  as  to  approximate  the  actual  unknown  current  distribution  reasonably 
well,  the  necessary  size  of  the  matrix  can  be  held  small  for  a  given  accu¬ 
racy  of  the  solution,  resulting  in  numerical  efficiency.  Moreover,  the 
numbers  of  basis  functions  for  y/z)  and  JJiz)  should  be  equal  [Ref.  18]. 
Hence  the  choice  of  basis  functions  is  important  from  the  numerical  point 
of  view.  In  actual  computations  for  the  dominant  mode  [Ref.  1],  J^,  and 
7,,  have  been  chosen  to  be  product  functions  of  the  form: 

(6a) 
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ut  m 


Note  that  (6a, b)  are  the  Fourier  transforms  of 
J2iix,z)  =  Ji{x)J2(z) 

Jxl(x,z)  =  J^(x)J4(z) 

Where  the  functional  forms  of  J\  y  Ji  >  ,  Ji  ,  are  given  in  Figure  2. 

The  coordinate  forms  of  the  current  distributions  shown  in  Figure  2  are 
the  following: 


J2{z)  =  y  COS(  -y ) 

(76) 

sin(  ^ ) 

(7c) 

- 

[Id) 
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Taking  the  Fourier  transform  of  (7)  we  get  : 


J\K)  = 


2sin(ct^w) 


ct„w 


+ 


[  cos(a„w)  - 


2vin(a„w)  2{1  -  cos(a„w)} 


+ 


]  (8a) 


= 


TT  cos(j3/) 


(m 
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Figure  2.  Forms  of  assumed  current  distributions. 


^3(0 


In  sin(a„iv) 


(8c) 


yAiu = 


cos{(iP} 

~ir 


sin(^/) 


m 


Corresponding  to  the  structure  in  Figure  1,  w  in  (7)  and  (8)  is  the  half 
width  of  the  strip  and  /  is  the  half  length  of  the  strip. 
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2.  Procedure  of  solution  for  resonant  frequency 
The  field  equations  (la,b) 

^zr^z  ^zx^x  ~  ^z 

were  derived  assuming  the  existence  of  ,  and  any  solution  derived  from 
(la,b)  must  satisfy  the  equations  (la)  and  (lb)  simultaneously.  Since  there 
are  four  unknowns,  Cj  ,  ,  Ex  in  (la,b),  they  can  not  be  solved  in 

general.  We  introduce  Galerkin's  procedure  by  taking  inner  products  of 
(la)  with  J^i  and  likewise  (lb)  with  Jxt  .  Then,  apply  the  (3a,b), 


By  the  virtue  of  Parseval's  theorem,  the  right  sides  of  (9)  are  zero,  so  we 
have 


(10a) 
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In  second-order  approximation,  equations  (10)  form  the  matrix  solution 
with  dimensions  [2x2]  .  The  matrix  is  symmetrical.  In  order  to  have 

nontrivial  values  for  Cj  ,  d,^  ,  the  determinant  of  the  matrix  must  be  zero. 
a.  First  order  approximation 

For  the  lowest-order  approximation,  here  called  first-order,  we 
assume  the  current  only  in  z-direction  ^  =  UJ  =  1),  and  choose  the 

current  distribution  =  y,(x)  x  /jW-  Since  we  have  current  only  in  the 
z-direction  (  Jx  =  0  ),  then  from  (10)  there  remains  only  the 

equation: 


^^1 


Z 


/z1'^2Z‘4i 


=  0 


a,P 


(11) 


By  using  trial  values  of  frequency,  the  solution  (  resonant  frequency  ) 

is  found  as  the  value  best  satisfying  (11) 
b.  Second  order  approximation 

For  the  second  order  solution  we  assume  current  exists  in  both 
directions  (  z  and  x  )  and  choose  the  current  distribution. 


/rl  (®/i,^m)  “ 
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Since  we  have  current  in  the  z  and  x-directions  (  7^.  #  0,  0),  then  we 

have  the  matrix  solution 


<^1  7  ^zX^ir^zX  +  7  ^zX^ix^xX  “  ^ 


a,  ^ 


a,  P 


(12a) 


Cl 


Txl'^xzTzl 


+ 


a,/? 


a,P 


(126) 


In  this  case  we  find^  as  the  frequency  for  which  the  determinant  of  the 
coefficients  vanishes. 

C  METHOD  OF  SOLUTION  FOR  PROPAGATION  CONSTANT 

The  increasing  use  of  microstrip  lines  at  microwave  frequencies  has 
recently  created  considerable  interest  in  the  study  of  dispersion  charac¬ 
teristics  of  these  lines.  Until  quite  recently,  much  of  the  work  on  the 
microstrip  line  was  based  on  a  TEM  analysis.  This  analysis  is  employed 
to  calculate  the  static  capacitance  of  the  structure  from  which  the  char¬ 
acteristic  impedance  and  the  propagation  wavenumber  are  subsequently 
derived.  However,  this  analysis,  which  is  necessarily  approximate,  is  inef¬ 
fective  for  estimating  the  dispersion  properties  of  the  line  at  higher  fre¬ 
quencies  [Ref.  14]. 

1.  Field  equation 

The  analysis  technique  for  computing  dispersion  characteristics  of 
shielded  transmission  line  of  the  type  considered  here  has  been  carried 
out  by  Mittra  and  itoh  [Ref.  3].  The  formulation  of  the  problem  applied 
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and  Jj.  are  the  transforms  of  strip  currents  J,  and  J^.  Notice  that  and 
Ej.  are  unknown  since  the  electric  fields  E,{x,D  +  T)  and  E^{x,D  +  T)  are 
unknown  for  vv  <  |  x  |  <  a,  though  thev  are  zero  on  the  strip.  The  field 
equations  in  (13a,b)  are  similar  to  the  field  equations  in  (la,b),  but  with 
current  distributions  assumed  uniform  in  the  z-direction  and  propagation 
constant  P  is  now  the  unknown. 
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2.  Procedure  of  solution  for  propagation  constant 

The  first  step  is  to  expand  the  unknowns  and  in  terms  of 
known  basis  functions  J^j  and  7^^* 

=  (15^) 

J=1 

N\ 

^=1 

Where  Cj  and  ^4  ^^e  unknown  coefficients.  The  basis  functions  and 
J^j  must  be  chosen  such  that  their  inverse  Fourier  transforms  are  nonzero 
only  on  the  strip  |  x  j  <  w.  After  substituting  (15)  into  (13)  and  taking 
inner  products  with  the  basis  functions  J^t  and  7^,  for  different  values  of 
i.  This  process  yields  the  equation, 


A'l 

A'2 

II 

(16u) 

7=1 

k=l 

A'l 

N2 

=  0 

,1=  l,2,3,..,iVl. 

(165) 

7=1 

A=1 

Where  from  definition  of  inner  products  associated  with  the  Fourier 
transform  defined  by  (14),  the  matrix  elements  are 


17 


4''"®  =  9  ^  ( I 

n=—oo 

+  00 

4‘'^*OT  =  rf*  ^  JMZ.A'^nPJjxkK)  (>7« 

«=— oo 

+  00 

4^’'\/J)  =  9  ^  7„(«„)Z„(«„;„,)7,/or„)  (17c) 

M=  — OO 

+CO 

«=  — OO 

The  right  hand  side  of  (13)  can  be  eliminated  in  Galerkin's  process  via 
the  application  of  Parseval's  relation. 

Now  the  simultaneous  equations  (16)  are  solved  for  the  propa¬ 
gation  constant  P  at  each  frequency  co  by  setting  the  determinant  of  the 
coefficient  matrix  equal  to  zero. 

a.  First  order  approximation. 

We  choose  of  matrix  size  ;  N1  =0,  N2=  1.  The  matrix  problem 
reduces  to  ; 


c,yi,(a)4(“.ffi^HW  =  0  (18) 
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b.  Second  order  approximation. 

We  choose  matrix  of  size  :  N1  =  1,  N2=  1.  The  matrix  to  be 

solved  is 

C|yy,|(a)Z,,(c<,  =  0  (19a) 

a,/?  a,p 

c,y  7,i(a)Z„(<ic,  P)7,,(oc)  +  </,y  7^,(a)Z,^{«,  P)7^,(«)  =  0  (19*) 

a,/?  a,p 

Typical  output  values  for  the  moment-method  calculation  of 
propagation  constant  P  for  various  frequencies  are  shown  in  Table  3  and 
Table  4  on  page  24,  Figure  6  on  page  29  and  Figure  7  on  page  30. 

D.  METHOD  OF  SOLUTION  FOR  FRINGING  CAPACITANCE 

In  this  work,  the  resonant  frequency  data  of  the  suspended  substrate 
resonators  have  been  used  for  estimating  the  end  effect  at  the  open  ends 
of  suspended  substrate  lines  using  the  full-wave  theory.  In  order  to  get  the 
value  of  fringing  capacitance,  we  use  two  methods,  the  equivalent  length 
extension  model  and  the  transmission  line  circuit  model. 

The  descriptions  of  each  method  are  as  follows  : 

1.  Equivalent  length  extension  model 

From  the  dispersion  relation  the  guide  wavelength  Xg  =  2-^  is  de¬ 
rived  at  the  resonant  frequency  ^  of  the  microstrip  resonator  of  length  I. 
Consider  the  open-circuited  resonator  whose  structure  is  identical  to 
Figure  1  except  that  the  length  is  /'  instead  of  /.  The  length  /'  is  deter- 
mined  from  the  resonant  condition  of  the  open  circuited  line  /'  =  — . 
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The  hypothetical  extension  of  the  suspended  substrate  line,  which  ac¬ 
counts  for  the  end  effect,  is  given  by  A/  =  /'  —  /.  In  terms  of  TEM  model 
for  transmission  line  characteristic  impedance  is  defined  by  L 
(Henry/meter)  and  Q  (Farad/ meter),  and  we  have  equations; 

z„  =  vi/Q  W 

p  =  (20b) 


From  equations  (20a)  and  (20b),  the  capacitance  C„  is  defined 


(21) 


where  is  capacitance/meter  of  line.  The  value  of  fringing  capacitance 
AC  is 


AC  =  A/Co  =  (/'-/) 


P 


In  this  work  we  calculate  the  end  effect  for  a  single  open  end,  as  shown 
in  Figure  3  (  for  )  using  the  computer  program  in  Appendix  D. 

2.  Transmission  line  circuit  model 

The  value  of  fringing  capacitance  can  be  calculated  by  using  the 
approach  of  the  transmission  circuit  model.  The  input  impedance  seen 
from  the  center  of  the  strip  resonator  is  calculated  using  formula  as  the 
following  [Ref.  12]: 


Zc+jZ^idnipt) 
Z^  +JZ,  taniPD 


(23) 
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MlCKOirrJII^  X£5O^a70> 

1 

4 

1 

Figure  3.  End  effect  on  substrate  resonator 
Where 

2  s= _ 1 _ 

ycu^AC 

and  AC  is  a  fringing  capacitance,  /  is  the  half  of  strip  length,  P  is  the 
propagation  constant,  obtained  as  explain  in  the  previous  section,  and 
Zo  is  the  characteristic  impedance  of  strip.  The  formula  to  be  used  to 
calculate  characteristic  impedance  Z„  is  the  following  [Ref.  17]: 

For  Q<W<a 


Zi  =  60[  F  +  R  /„{6G/  W  +  J\+4iGIW)^  }] 
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Where, 


V=  -1.7866  -  0.2035(r/G)  +  0.4750(2^/(7) 

R  =  1.0835  +  0.1007(7/0  -  0.09457(2/1/0 

/ - _ 

y^eff  —  cj 

and  c  is  the  free-space  velocity  of  light.  This  formula  has  been  confirmed 
by  variational  calculation  [Ref.  19],  and  is  assumed  here  to  vary  only 
slowly  with  frequency.  The  computer  program  for  this  calculation  is 
shown  in  Appendix  D.  In  Table  2  is  shown  the  values  of  Z^,  the  char¬ 
acteristic  impedance  for  different  value  of  and  the  width  of  strip  . 


Table  2.  THE  CHARACTERISTIC  IMPEDANCE  OF  SUSPENDED 
SUBSTRATE  LINE 


Dielectric 

constant 

(£r) 

Width  of 

strip 

(mm) 

Impedance  first  order 
(jS).  Zo 

Impedance  second 
order  (j5), 

2.2 

0.711 

94.635 

94.220 

1.0 

79.185 

78.789 

1.25 

68.963 

68.596 

9.8 

0.711 

70.651 

70.204 

1.0 

60.533 

59.916 

1.25 

53.446 

52.732 

The  true  value  of  fringing  capacitance  AC  can  be  founded  by  using 
trial  values  in  (23)  and  looking  for  the  minimum  value  of  the  input 
impedance  Zj„. 
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III.  DISCUSSION  OF  RESONATOR  COMPUTATION 


A.  INTRODUCTION 

The  calculation  of  the  suspended  stripline  resonator  in  this  work  used 
numerical  computations  in  a  scale  model  of  millimeter  wave  integrated 
circuits.  The  numerical  results  presented  in  this  paper  have  been  gener¬ 
ated  on  a  main  frame  IBM  370  using  FORTRAN  77  and  compiler  VS  2 
FORTRAN,  since  this  compiler  is  10  times  faster  than  WATFOR  77. 
The  solution  for  resonant  frequencies  and  propagation  constant  was  made 
by  taking  the  absolute  value  of  the  determinant  matrix.  This  approach 
leads  to  the  minimum  value  for  the  system  determinant  rather  than 
crossing  the  axis  of  determinant  =  0.  All  output  data  together  are  shov.n 
in  Appendix  A,  while  the  complete  data  values  are  shown  in  Table  3  for 
£,  =  2.2  and  Table  4  for  e,  =  9.8.  The  computation  results  for  fringing 
capacitance  using  the  transmission  line  circuit  model  agree  with  the 
computation  result  using  the  length  extension  model,  as  is  shown  by  the 
very  small  difference  between  results  of  these  two  methods. 

For  each  of  these  methods  we  can  see  that  as  the  width  of  strip  in¬ 
creases  the  resonant  frequency  of  a  resonator  of  given  length  increases  as 
shown  in  Figure  4  and  Figure  5.  The  complex  dispersion  behavior  of  this 
line  system  is  shown  by  the  decrease  of  with  increasing  the  width  of 
strip.  Also  the  value  of  fringing  capacitance  goes  up  with  width.  This 
situation  may  be  due  to  the  fact  that  the  end  region  will  store  more  ac¬ 
cumulated  charge  with  wider  lines  of  the  strip.  There  will  be  currents 
flowing  in  the  end  region,  corresponding  to  the  extra  charge  [Ref.  6]. 
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Table  3.  THE  DATA  RESULTS  FOR  «,  =  2.2 


Method 
of  sol¬ 
ution 

Width 

of 

strip, 

w 

(mm) 

Reso¬ 

nant 

fre¬ 

quency, 

/o 

(GHZ) 

Propa¬ 
gation 
con¬ 
stant,  [i 
(Rad/m) 

End- 

effect 

model, 

AC, 

(Farad) 

Trans¬ 

mission 

line 

model, 

AQ 

(Farad) 

Error, 

AC„ 

First 

order 

0.711 

12.9096 

301.8001 

0.8050E-14 

0.8060E-14 

l.E-17 

1.0 

13.0237 

301.0026 

1.0152E-14 

1.0167E-14 

l.E-17 

1.25 

13.0913 

300.5175 

1.2024E-14 

1.2043E-14 

2.E-17 

■ 

0.711 

12.8582 

301.9217 

0.8038E-14 

0.8048E-14 

l.E-17 

1.0 

12.9668 

301.1930 

1.0099E-14 

1.0114E-14 

l.E-17 

1.25 

13.0331 

300.7837 

1.1906E-14 

1.1923E-14 

2.E-17 

Table  4.  THE  DATA  RESULTS  FOR  =  9.8 


Method 
of  sol¬ 
ution 

Width 

of 

strip, 

w 

(mm) 

Reso¬ 

nant 

fre¬ 

quency, 

/o 

(GHZ) 

Propa¬ 
gation 
con¬ 
stant,  (} 
(Rad/m) 

Trans¬ 

mission 

line 

model, 

ACj 

(Farad) 

Error, 

AC„ 

First 

order 

0.711 

9.4069 

294.5652 

2.3461E-14 

2.3537E-14 

8.E-17 

1.0 

9.7084 

293.5173 

2.795  IE- 14 

2.8050E-14 

lOE-17 

1.25 

9.8917 

292.9945 

3.1858E-14 

3.1977E-14 

l.E-16 

Second 

order 

0.711 

9.35427 

294.7848 

2.3477E-14 

2.3551E-14 

7.E-17 

1.0 

9.6251 

293.9983 

2.7820E-14 

2.7915E-14 

9.E-17 

1.25 

9.7804 

293.6181 

3.1695E-14 

3.1806E-14 

l.E-16 
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B.  RESONANT  FREQUENCY 

For  the  first  order  approximation  the  modified  summation  i'h*  terms 
of  a  and  to  (1 1)  is  represented  as 

+00  +00 

^  Y,  G(«„,W  =  0  (24) 

n=— oom=— oo 

Where  G(a„,  is  the  product  form  of  basis  functions  and  Green's 
impedance  function.  The  computer  programs  for  first  and  second  order 
approximation  for  is  shown  in  Appendix  B.  These  programs  calculate 
the  value  of  the  determinant  with  various  values  of  strip  width  in  order 
to  get  the  resonant  frequency,  with  dielectric  constant  for  substrate  : 

=  2.2  for  duroid,  and  =  9.8  for  alumina.  The  solutions  for  first  order 
are  found  : 

+30  +3000 

Y  X  ^l(“».  « =  0  (25) 

'’=-^^m=-3000 

In  the  x-direction  we  take  the  sum  from  from  n  =  -30  to  n  =  +30  and  in 
the  z-direction  from  m  =  -3000  to  m  =  +  3000.  Since  a  higher  limit  does 
not  change  the  result  as  shown  in  Table  5. 

In  a  similar  way,  we  may  compute  the  second  order  approximation 
with  the  same  summation  limits  for  a„  and  for  All  data  outputs  for 
resonant  frequency  appear  in  Appendix  A. 
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RESONANT  FREQUENCY 


THE  WIDTH  OF  STRIP  (MM) 


Figure  4.  Resonant  frequency  versus  width  of  strip  for  e,  -  2.2 


Table  5.  THE  COMPARISON  FOR  RESONANT  FREQUENCY  IN  FIRST 
ORDER  APPROXIMATION 


0-n 

w 

Resonant  Frequency 
(GHZ) 

n  =  -30  to 
n=  +30 

m  =  -3000  to 
m=  +3000 

2.2 

l.O 

13.023998 

n  =  -40  to 
n  =  +  40 

m  =  -4000  to 
m=  +4000 

2.2 

1.0 

13.023498 

RESONANT  FREQUENCY 


o 


I”!  _ 


0.70  0.76  0.82  0.88  0.94  1.00  1.06  1.12  1.18  1.24  1.30 


THE  WIDTH  OF  STRIP  (MM) 

Figure  5.  Resonant  frequency  versus  width  of  strip  for  s,  -  9.8 
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C.  PROPAGATION  CONSTANT 


The  program  to  be  used  for  computation  of  the  propagation  constant 
is  shown  in  Appendix  C.  An  initial  approximation  is  taken  as  the  fol¬ 
lowing  [Ref.  7]; 


(26) 


Where  c  is  the  free-space  velocity  of  light.  To  find  the  exact  value  of 
propagation  constant,  we  sweep  the  value  of  j5  from  1%  to  200%  of  (26) 
in  order  to  bring  the  value  of  the  matrix  determinant  in  (18)  to  zero.  The 
propagation  constant  can  be  identified  by  the  occurrence  of  a  minimum 
value  of  the  matrix  determinant.  In  this  case  we  take  the  summation 
limits  from  n  =  -40  to  n=  -^40,  since  a  higher  limit  does  not  change  the 
result  as  shown  in  Table  6. 


Table  6.  THE  COMPARISON  FOR  PROPAGATION  CONSTANT  IN  FIRST 
ORDER  APPROXIMATION 


w 

Propagation  Constant  ( 

n 

n  =  -40  to  n  =  -1-  40 

2.2 

1.0 

301.80 

o 

o 

+ 

II 

c 

o 

o 

o 

1 

II 

c 

2.2 

1.0 

301.81 

(  See  in  Appendix  A  for  the  complete  results  ) 
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PROPAGATION  CONSTANT 

300.0300.4  300.8  301.2  301.6  302.0 


FREQUENCY  (GHZ) 


Figure  6.  Resonant  frequency  versus  propagation  constant  for  e,  -  2.2 
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PROPAGATION  CONSTANT 

292.(202.98  293.  46  293.94  294.42  294.90 


FREQUENCY  (GHZ) 


Figure  7.  Resonant  frequency  versus  propagation  constant  fo”  e,  =  9.8 


D.  PREDICTION  OF  FRINGING  CAPACITANCE  VALUE 

The  value  of  fringing  capacitance  can  be  found  by  solving  equation 
(22)  for  the  equivalent  length  extension  model  and  solving  equation  (23) 
for  for  transmission  line  circuit  model.  The  computer  program  used  to 
calculate  the  fringing  capacitance  are  shown  in  Appendix  D.  From  the 
results  in  Table  3  and  Table  4  we  see  that  these  two  methods  agree  since 
ihere  is  only  a  small  difference  in  their  results. 


31 


FRINGING  CAPACITANCE  n0‘ 

8.5  9.5  10.5  11.5  12.5 


0.70  0.76  0.82  0.88  0.94  1.00  1.06  1.12  1.18  1.24  1.30 


THE  WIDTH  OF  STRIP  (MM) 


Figure  8.  Fringing  capacitance  versus  width  of  strip  for  £,  =  2.2 
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— I - 1 - 1 - 1 - 1  I  I  I 

0.82  0.88  0.94  1.00  1.06  1.12  1.18  1.24 

THE  WIDTH  OF  STRIP  (KIM) 


i  - 1 - 

0.70  0.76 


1.30 


Figure  9.  Fringing  capacitance  versus  width  of  strip  for  s,  =  9.8 


IV.  CONCLUSION 


From  a  computational  point  of  view,  the  calculation  using  Green's 
function  is  very  effective  to  solve  the  stripline  resonator  problem.  Actu¬ 
ally  the  computed  result  for  the  resonant  frequency  depends  on  the  choice 
of  the  definition  of  the  current  distributions  in  the  x-direction  and  the 
z-direction.  In  the  Table  3,  Table  4,  Figure  5  and  Figure  6  it  is  seen  that 
the  resonant  frequencies  increases  as  the  width  of  strip  increases.  We  do 
not  see  this  phenomenon  in  a  resonator  without  suspended  substrate,  as 
shown  in  [Ref.  1]. 

The  values  of  the  fringing  capacitance  found  here  appear  to  be  rea¬ 
sonable;  the  values  are  of  the  order  of  lO'"*  Farad.  The  error  between 
the  two  methods  of  calculation  is  very  small,  providing  a  confirmation 
of  the  length  extension  model. 
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APPENDIX  A.  OUTPUTDATA 


A.  RESONANT  FREQUENCY 
1.  First  order  approximation 
■  For  name  w  =  0. 711  mm  and  *,  =  2.  2 


FREQUENCY  (GHZ) 

12. 909532500000 
12. 909547800000 
12. 909563100000 
12. 909578300000 
12. 909593600000 


■  For  w  =  1. 0  mm  and 

FREQUENCY  (GHZ) 

13. 023622500000 
13. 023684500000 
13. 023746500000 
13. 023808500000 
13. 023870500000 


■  For  w  =  1.  25  mm  and 
FREQUENCY  (GHZ) 

13. 090996700000 
13. 091095900000 
13. 091195100000 
13. 091294300000 
13. 091393500000 
13. 091492700000 
13. 091591800000 


•  For  w  -  0. 711  mm  and  t 
FREQUENCY  (GHZ) 

9.406851964275 
9. 406858630475 
9.406865296675 
9.406871962875 
9. 406878629075 


DETERMINANT'S  VALUES 

0.0195708908 
0. 0134981312 
0. 0063938200 
0. 0007518493 
0.0081374571 


2.  2 

DETERMINANT’S  VALUES 

0.0549158342 
0. 0306886621 
0. 0076220483 
0.0151751861 
0. 0382801555 


=  2.2 

DETERMINANT'S  VALUES 

0. 0922133327 
0. 0595099851 
0. 0294705480 
0. 0022094622 
0. 0342780240 
0.0658574104 
0.0986036658 


=  9.8 

DETERMINANT'S  VALUES 

0.0055711373 
0. 0026056429 
0. 0003598495 
0.0033253400 
0.0062908285 
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•  For  w  =  1.  0  tnm  and  = 
FREQUENCY  (GHZ) 


9. 708079097500 
9. 708178995000 
9. 708278892500 
9. 708378790000 
9. 708478687500 
9. 708578585000 
9. 708678482500 


■  For  w  =  1.  25  nun  and 

FREQUENCY  (GHZ) 


9. 891719662500 
9. 891720994500 
9. 891722326500 
9. 891723658500 
9. 891724990500 
9. 891726322500 
9. 891727654500 


.  8 

DETERMINANT’S  VALUES 


0. 1101123587 
0.0734216453 
0. 0367312736 
0. 0000412434 
0. 0366484453 
0. 0733377923 
0. 1100267979 


9.8 

DETERMINANT'S  VALUES 


0.0013436906 
0. 0009188024 
0. 0004939143 
0. 0000690263 
0.0003558617 
0. 0007807497 
0.0012056375 


2.  Second  order  approximation 

■  For  N  =  0.711  mm  and  =  2.  2 

FREQUENCY  (GHZ)  THE  VALUES  OF  DETERMINANT 


12.  858237475000 
12.  858239430000 
12.  858241385000 
12.  858243340000 


134.  7679778482 
45. 1968072590 
44.3743231685 
133.9454133199 


■  For  w  =  1.  0  mm  and  =  2.  2 
FREQUENCY  (GHZ) 


12. 966809460000 
12. 966809956000 
12. 966810452000 
12.  966810948000 
12.  966811444000 


THE  VALUES  OF  DETERMINANT 


21.9744171435 
12. 5409965103 
3. 1075768752 
6. 3258416755 
15. 7592591259 


■  For  w  =  1.  25  mm  and  =  2.  2 

FREQUENCY(GHZ)  THE  VALUES  OF  DETERMINANT 


13.  033089704000 
13. 033102202000 
13.  033114700000 


95.  3389230655 
36.4297524373 
168.  1980546933 
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■  For 


■  For 


■  For 


=  0. 711  mm  and  =  9. 
FRZQUENCY(GHZ) 


9.  354219440000 
9. 354244230000 
9. 354269030000 
9.  354293820000 


=  1.  0  mm  and  =  9.  8 
FREQUENCY(GHZ) 


9. 624875000000 
9. 625000000000 
9. 625125000000 
9. 625250000000 
9. 625375000000 


=  1.  25  mm  and  =  9.  8 
FREQUENCY  (GHZ) 


9. 780300000000 
9. 780360000000 
9. 780420000000 
9. 780480000000 


THE  VALUES  OF  DETERMINANT 


838. 0000000000 
369. 8750000000 
96. 0000000000 
605. 0000000000 


THE  VALUES  OF  DETERMINANT 


2076. 6096794134 
1067. 1692908816 
57. 7672394809 
951.5964768065 
1960. 9218599920 


THE  VALUES  OF  DETERMINANT 


359. 7268520610 
83. 7447038466 
192. 2324952445 
468. 2047453591 
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B.  PROPAGATION  CONSTANT 
1.  First  order  approximation 


■  For  w  =  0.  711  mm  and 
FREQUENCY  (GHZ) 

=  2.2 

BETA  (RAD/METER) 

DETERMINANT 

12. 909578300 

301. 793213000 

0. 006551314 

12. 909578300 

301. 796631000 

0. 003314178 

12. 909578300 

301. 800049000 

0.  000016721 

12. 909578300 

301. 803467000 

0. 003176749 

12. 909578300 

301. 806885000 

0.  006300699 

■  For  !¥  —  1.  0  nun  and  = 

2.  2 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

13. 023746500 

300. 933559227 

0. 043851649 

13. 023746500 

300. 968062043 

0. 016535467 

13. 023746500 

301. 002564859 

0. 010783630 

13.  023746500 

301.  037067675 

0.  038105642 

13. 023746500 

301. 071570490 

0.  065430568 

■  For  w  =  1. 25  mm  and  = 

=  2.2 

IREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

13. 091294300 

300.448130204 

0. 047616359 

13. 091294300 

300. 482811969 

0. 023635390 

13. 091294300 

300.517493734 

0. 000348146 

13. 091294300 

300. 552175499 

0. 024334248 

13. 091294300 

300.  586857264 

0. 048322918 

■  For  w  —  0.  711  mm  and 

=  9.8 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

9. 406865000 

294. 290520504 

0. 218513360 

9. 406865000 

294. 427868289 

0.  122327343 

9. 406865000 

294. 565216074 

0. 026102604 

9. 406865000 

294. 702563859 

0. 070160849 

9. 406865000 

294. 839911643 

0.  166463009 

•  For  w  =  1.  0  mm  and  = 

9.  8 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

9. 708379000 

293. 233772058 

0.  208247936 

9. 708379000 

293. 375522190 

0. 123411162 

9. 708379000 

293. 517272321 

0. 038539359 

9. 708379000 

293. 659022453 

0. 046367465 

9. 708379000 

293. 800772584 

0. 131309303 
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9. 708379000 


293. 942522716 


0. 216286146 


For 

w  =  1.  25  mm  and  = 

FREQUENCY  (GHZ) 

9.  8 

BETA  (RAD/METER) 

DETERMINANT 

9. 891724000 

292. 705621320 

0. 146004764 

9. 891724000 

292. 850048436 

0. 069357733 

9. 891724000 

292. 994475552 

0. 007321405 

9. 891724000 

293. 138902668 

0. 084032642 

9, 891724000 

293. 283329783 

0.  160775971 

2. 

Second  order  approximation 

For 

w  =  0. 711  mm  and 

=  2.2 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

12. 858241000 

301. 914079040 

264. 8257310433 

12. 858241000 

301. 916633866 

188. 2549021972 

12. 858241000 

301. 919188693 

111.  6834688237 

12. 858241000 

301. 921743520 

35. 1114309231 

12. 858241000 

301. 924298346 

41.  4612114998 

12. 858241000 

301. 926853173 

118. 0344584496 

12. 858241000 

301. 929407999 

194. 6083099229 

12. 858241000 

301.931962826 

271. 1827659171 

For  w  =  1.  0  mm  and  = 

2.  2 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

12. 966810000 

301. 189557419 

47. 3499570198 

12. 966810000 

301. 191275018 

25. 5508907438 

12. 966810000 

301. 192992616 

3. 7517095664 

12. 966810000 

301. 194710215 

18. 0475865116 

12. 966810000 

301. 196427814 

39. 8469974898 

For  w  -  1.  25  mm  and  tj. 

=  2.2 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

13. 033102000 

300.  776834261 

27. 3862465221 

13. 033102000 

300. 778560641 

15. 0664900067 

13. 033102000 

300. 780287021 

2. 7466685987 

13. 033102000 

300. 782013401 

9. 5732177036 

13. 033102000 

300. 783739781 

21. 8931688989 

For  w  =  0. 711  mm  and 

=  9.8 

FREQUENCY  (GHZ) 

BETA  (RAD/METER) 

DETERMINANT 

9. 354269000 

294. 739297051 

402.  7507924038 

9. 354269000 

294.  762060358 

192. 4514179721 
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9.354269000  294.784823665 
9.354269000  294.807586971 
9.354269000  294.830350278 


■  For  w  =  1.  0  mm  and  =  9.8 

FREQUENCY  (GHZ)  BETA  (RAD/METER) 


9.625125000  293.904601953 
9.625125000  293.951446805 
9.625125000  293.998291657 
9.625125000  294.045136509 
9.625125000  294.091981361 


■  For  w  —  1.  25  mm  and 

FREQUENCY  (GHZ) 


=  9.  8 

BETA  (RAD/METER) 


9. 780360000 
9. 780360000 
9. 780360000 
9. 780360000 
9. 780360000 
9. 780360000 
9. 780360000 


293. 532445612 
293.  561005835 
293. 589566057 
293. 618126279 
293  646686502 
293. 675246724 
293. 703806947 


17. 8618972464 
228. 1891527630 
438. 5303480889 


DETERMINANT 


311. 5213476632 
120. 8819158085 
69. 7832174509 
260. 4740500634 
451. 1905799777 


DETERMINANT 


204. 6350496775 
136. 5694649559 
68. 4983401777 
0. 4216756325 
67. 6605283907 
135. 7482716025 
203. 8415537136 
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C.  FRINGING  CAPACITANCE  USING  LENGTH  EXTENSION  MODEL 

1 .  First  order  approximation 

■  For  w  =  0.  711  mm  and  Cjr  —  2.  2 

RESONANT  PROPAGATION  FRINGING*'CAPACITANCE 

FREQUENCY  CONSTANT  (  FARAAD  ) 

12.9095780  301.800050  0.00000000000000805039 


•  For  AT  =  1.  0  mm  and  =  2.  2 

RESONANT  PROPAGATION  FRINGING  CAPACITANCE 

FREQUENCY  CONSTANT  (  FARAAD  ) 

13.0237470  301.002570  0.00000000000001015215 


•  For  w  =  1.  25  mm  and  =  2.  2 

RESONANT  PROPAGATION  FRINGING  CAPACITANCE 

FREQUENCY  CONSTANT  (  FARAAD  ) 

13.0912940  300.517490  0.0000000000000120243 


■  For  w  =  0. 711  mm  and  =  9.  8 

RESONANT  PROPAGATION  FRINGING  CAPACITANCE 

FREQUENCY  CONSTANT  (  FARAAD  ) 

9.4068650  294.565220  0.000000000000023461 


■  For  =  1.  0  mm  and  =  9.  8 

RESONANT  PROPAGATION  FRINGING  CAPACITANCE 

FREQUENCY  CONSTANT  (  FARAAD  ) 

9.7083790  293.517270  0.000000000000027951 


■  For  w  =  1.  25  mm  and  =  9.8 

RESONANT  PROPAGATION  FRINGING  CAPACITANCE 

FREQUENCY  CONSTANT  (  FARAAD  ) 

9.8917240  292.994480  0.00000000000003185787 
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2.  Second  order  approximation 

■  For  w  —  0.  711  nun  and  =  2.  2 

RESONANT  PROPAGATION 

FREQUENCY  CONSTANT 


12.8582410  301.921740 


FRINGING  CAPACITANCE 
(  FARAAD  ) 

0. 0000000000000080382 


For  V  =  1.  0  mm  and  =  2.  2 

RESONANT  PROPAGATION 

FREQUENCY  CONSTANT 


12. 9668100 


301. 192990 


FRINGING  CAPACITANCE 
(  FARAAD  ) 

0. 000000000000010099657 


For  w  =  1.  25  mm  and  tj-  =  2.  2 

RESONANT  PROPAGATION 

FREQUENCY  CONSTANT 


13. 0331020 


300. 783740 


For  w  =  0.711  mm  and  cj-  “  ® 

RESONANT  PROPAGATION 

FREQUENCY  CONSTANT 


9. 3542690 


294. 784800 


For  w  =  1.  0  mm  5ind  =  9.  8 

RESONANT  PROPAGATION 

FREQUENCY  CONSTANT 


9. 6251250 


293. 998290 


For  w  =  1.  25  mm  and  -  9.8 

RESONANT  PROPAGATION 

FREQUENCY  CONSTANT 


9. 7803600 


293. 618130 


FRINGING  CAPACITANCE 
(  FARAAD  ) 

0. 0000000000000119056853476236 


FRINGING  CAPACITANCE 
(  FARAAD  ) 


0. 0000000000000234773 


FRINGING  CAPACITANCE 
(  FARAAD  ) 

0. 0000000000000278197027432552 


FRINGING  CAPACITANCE 
(  FARAAD  ) 


0. 0000000000000316946 
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D.  FRINGING  CAPACITANCE  USING  TRANSMISSION  LINE  CIRCUIT 

• 

MODEL 

1.  First  order  approximation 

■  For  w  -  0. 711  mm  and 

=  2.2 

FREQUENCY  (GHZ) 

CAPACITANCE  (FARAAD) 

IhPEDANCE 

12. 9095783 

0. 00000000000000800 

0. 043897341449 

12. 9095783 

0. 00000000000000802 

0, 029423400665 

12. 9095783 

0. 00000000000000804 

0. 014949733818 

12. 909578’ 

0. 00000000000000806 

0. 000476340899 

12. 9095783 

0. 00000000000000808 

0. 013996778098 

12. 9095783 

0. 00000000000000810 

0. 028469623184 

12. 9095783 

0. 00000000000000812 

0. 042942194364 

• 

■  For  w  =  1.  0  mm  and  tj.  = 

2.  2 

FREQUENCY  (GHZ) 

CAPACITANCE  (FARAAD) 

IMPEDANCE 

13. 0237470 

0. 0000000000000101667 

0. 000020988378 

13. 0237470 

0. 0000000000000101668 

0. 000010770731 

13. 0237470 

0, 0000000000000101668 

0. 000000553083 

13. 0237470 

0. 0000000000000101668 

0. 000009664563 

13. 0237470 

0. 0000000000000101668 

0. 000019882210 

13.0237470 

0. 0000000000000101668 

0. 000030099857 

13.0237470 

0. 0000000000000101669 

0. 000040317503 

■  For  w  =  1.  25  mm  and  tj.  = 

=  2.2 

FREQUENCY  (GHZ) 

CAPACITANCE  (FARAAD) 

IMPEDANCE 

13, 0912940 

0.  000000000000012043 

0. 0000123459405 

13. 0912940 

0.  000000000000012043 

0. 0000045583584 

13. 0912940 

0.  000000000000012043 

0. 0000032292234 

13, 0912940 

0.  00000 jOOOOOOO 12043 

0. 0000110168053 

13. 0912940 

0.  000000000000012043 

0. 0000188043870 

•  For  w  =  0,  711  mm  and 

=  9.8 

FREQUENCY  (.GHZ) 

CAPACITANCE  (FARAAD) 

IMPEDANCE 

9.  4068650 

0.  00000000000002353658 

0. 000004112742 

9. 4068650 

0.  00000000000002353658 

0. 000002943926 

9.4068650 

0. 00000000000002353658 

0. 000001775109 

9. 4068650 

0.  00000000000002353659 

0. 000000606293 

9.4068650 

0. 00000000000002353659 

0. 000000562522 

9. 4068650 

0. 00000000000002353660 

0. 000001731339 

9. 4068650 

0. 00000000000002353660 

0. 000002900155 

■  For  nr  =  1.  0  mm  and  = 

9.  8 
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FREQUENCY  (GHZ)  CAPACITANCE  (FARAAD'' 


9. 7083790 
9. 7083790 
9. 7083790 
9. 7083790 
9. 7083790 
9. 7083790 


0. OOOOOOOOOOOOC28040 
0. 000000000000028045 
0. 000000000000028050 
0. 000000000000028055 
0. OOU000000000028060 
0. 000000000000028065 


■  For  w  =  1.  25  mm  and  tj- 
FREQUENCY  (GHZ) 


9. 8917240 
9. 8917240 
9. 8917240 
9. 8917240 
9. 8917240 


=  9.  8 

CAPACITANCE  (FARAAD) 


0. 000000000000031975 
0. 000000000000031976 
0. 000000000000031977 
0. 000000000000031978 
0. 000000000000031979 


2.  Second  order  approximation 


■  For  J/  =  0.  711  ram  and  —  1.2. 

FliEQUENCY  (GHZ)  CAPACITANCE  (FARAAD) 


12.  8582410 
12.  8582410 
12. 8582410 
12,  8582410 
12. 8582410 


0.  000000000000008044 
0. 000000000900008046 
0. 000000000000008048 
0. OOOOOOOOOOOCD08050 
0. 000000000000008052 


■  For  v  -  1,  0  mm  and  tj-  — 
FREQUENCY  (GHZ) 


12,  9668100 
12. 9668100 
12. 9668100 
12.  9668100 
12. 9668100 
12. 9668100 


2.  2 

CAPACITANCE  (FARAAD) 


0. 000000000000010111 
0. 000000000000010112 
0. 000000000000010113 
0.  0000000000v'»010114 
0. 000000000000010115 
0.  000000000000010116 


■  For  w  =  1.  25  mm  and  tj. 

FREQUENCY  (GHZ) 


13.  0331020 
13. 0331020 
13. 0331020 
13. 0331020 
13. 0331020 
13. 0331020 
13. 0331020 


=  2.2 

CAPACITANCE  (FARAAD) 


0. 00000000000001192344 
0.  00000000000001192345 
0. 00000000000001192346 
0.  00000000000001192347 
0. 00000000000001192348 
0. 00000000000001192349 
0. 00000900000091192350 


IMPEDANCE 


0. 002439262912 
0.  00133'’532326 
0. 000227305924 
0. 000877916293 
0. 001983634327 
0. 003089348177 


IMPEDANCE 


0. 000409506598 
0. 000233952816 
0. 000058399157 
0. 000117154378 
0. 000292707792 


IMPEDANCE 


0. 003037869907 
0. 001608804997 
0. 000179742742 
0. 001249316856 
0. 002678373799 


IMPEDANCE 


0. 001425888209 
0. 000922250353 
0. 000418612916 
0. 000085024103 
0. 000588660705 
0. 001092296888 


IMPEDANCE 


0. 000010349874 
0. 000006513849 
0. 000002677825 
0. 000001158198 
0. 000004994223 
0. 000008830247 
0. 000011666272 
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■  For 


■  For 


■  For 


=  0.  711  mm  and  Ej. 

=  9.  8 

FREQUENCY 

(GHZ) 

CAPACITANCE  (FARAAD) 

IMPEDANCE 

9. 3542690 

0.  0000000000000235503 

0. 000205384374 

9. 3542690 

0.  0000000000000235506 

0. 000119294322 

9. 3542690 

0. 0000000000000235509 

0. 000033204290 

9. 3542690 

0.  0000000000000235512 

0.  000052885720 

9. 3542690 

0.  0000000000000235515 

0.  000138975711 

=  1.  0  mm  and  tj.  = 

9.  8 

FREQUENCY 

(GHZ) 

CAPACITANCE  (FARAAD) 

IMPEDANCE 

9. 6251250 

0. 0000000000000279130 

0.  000283333299 

9. 6251250 

0. 0000000000000279135 

0.  000175879548 

9. 6251250 

0.  0000000000000279140 

0.  000068425835 

9. 6251250 

0. 0000000000000279145 

0.  000039027837 

9. 6251250 

0. 0000000000000279150 

0.  000146481471 

9. 6251250 

0. 0000000000000279155 

0.  000253935066 

9. 6251250 

0. 0000000000000279160 

0.  000361388622 

=  1.  25  mm 

and  Sj.  • 

=  9.8 

FREQUENCY 

(GHZ) 

CAPACITANCE  (FARAAD) 

IMPEDANCE 

9. 7803600 

0.  0000000000000318055 

0.  000174399228 

9. 7803600 

0.  0000000000000318058 

0.  000123674994 

9. 7803600 

0. 0000000000000318061 

0.  000072950771 

9. 7803600 

0. 0000000000000318064 

0.  000022226558 

9. 7803600 

0. 0000000000000318067 

0.  000028497644 

9. 7803600 

0. 0000000000000318070 

0.  000079221837 

9. 7803600 

0. 0000000000000318073 

0.  000129946020 
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APPENDIX  B.  PROGRAMS  FOR  CALCULATION  OF  RESONANT 

FREQUENCY 


•k 

* 
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* 

* 

* 

* 
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* 
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* 

* 

* 

* 

* 

* 

* 

* 

* 
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* 
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* 

* 

* 

* 

* 

* 

* 
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■A 

A 
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VARIABLE  DEFINITION: 


I  =  REPRESENT  THE  REGION  1,2,3  IN  RESONATOR 

ALPA  =  DISCRETE  TRANSFORM  VARIABLE  FOR  DOMINANT  MODE  (x-direction) 

BETA  =  DISCRETE  TRANSFORM  VARIABLE  FOR  DOMINANT  MODE  (z-direction) 

A  =  A  HALF  WIDTH  OF  RESONATOR  DEVICE  =  1.6  mm 

D  =  HEIGHT  OF  RESONATOR  BELOW  SUBSTRATE  =  0.  66  mm 

T  =  THE  THICKNESS  OF  SUBSTRATE  =  0. 255  mm 

W  =  THE  WIDTH  OF  STRIP,  CAN  BE  0.711  mm,  1.0  mm,  or  1.25  mm 

H  =  THE  HEIGHT  OF  RESONATOR  ABOVE  SUBSTRATE  =  0.  66  mm 

AL  =  THE  LENGTH  OF  STRIPLINE  ,  ASSUME  =  10  mm 

AKI  =  RADIAN  FREQUENCY  *  SQRT  OF  (PERMITTIVITY*PERMEABILITY) 

MO  =  PERMEABILITY  IN  VACCUM  =  PI*4.E-7 

MR  =  PERMEABILITY  IN  A  SUBSTRATE  ASSUME  =  1 

EO  =  PERMITTIVITY  IN  VACCUM  =  ( l/( 36*PI) )*E-9  =  8.85E-12 

ER  =  PERMITIVITY  IN  A  SUBSTRATE,  CAN  BE  2.2  OR  9.8 

FREQ  =  FREQUENCY  OPERATION  IN  GHZ 

DE  =  DENOMINATOR  OF  ZE 

DE  =  NUMERATOR  OF  ZE 

DH1,DH2,DH  =  DENOMINATOR  OF  ZH 

NH  =  NUMERATOR  OF  ZH 

ZE  =  IMPEDANCE  OF  ELECTRIC  FIELD,  eqn  ( Ig) 

ZH  =  IMPEDANCE  OF  MAGNETIC  FIELD,  eqn  (Ih) 

ZZZ,ZZX,ZXZ,AND  ZXX  ARE  THE  GREEN'S  IMPEDANCE  FUNCTION  (lc,d,e,f) 
COTl  =  CTl,  EQN  (li) 

COT2  =  CT2,  EQN  (Ij) 

C0T3  =  CT3,  EQN  (Ik) 

J1  =  EQN  (8a) 

J2  =  EQN  (8b) 

JJl  =  Jzl  EQN  (6a) 

JJ4  =  Jxl  EQN  (6b) 

JZl  =  BASIS  FUNCTION  OF  CURRENT  COMPONENT  IN  Z  -  DIRECTION 

JXl  =  BASIS  FUNCTION  OF  CURRENT  COMPONENT  IN  X  >  DIRECTION 

K11,K12,K21,K22  =  MATRIX  ELEMENT 

Jl,  J3  =  ASSUME  CURRENT  DISTRIBUTION  IN  X  -  DIRECTION 
J2,  J4  =  ASSUME  CURRENT  DISTRIBUTION  IN  Z  -  DIRECTION 
ALL  COMPUTATION  IN  STANDARD  "MKS" 
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A.  FIRST  ORDER  APPROXIMATION 

*  THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  THE  VALUE 

*  OF  DETERMINANT  OF  GREEN  FUNCTION  IN  FIRST  ORDER 

*  APPROXIMATION  WITH  W  =  0.711  mm  and  ER  =  2.  2 

*  .  . . . .  .  -  -  -■  ,  ■ 


*  MAIN  PROGRAM  : 

*  — 


*  FUNCTION  :  Fl( ALPA.BETA.FREQ)  =  ZZZ  *  JZl  **2 
COMPLEX  KK,K11,F1 

INTEGER  N,M,L,I 

REAL  A,W,D,T,H,MO,EO,ER,PI,MAGKK,AL,FREQ,ALPA,BETA 
PARAMETER  (PI  =  3.141592654,  EO  =  8.85E-12) 

PARAMETER  (H  =  0. 66E-03,  T  =  0. 255E-03,  D  =  0. 66E-03) 
PARAMETER  (MO  =  PI*4.E-7,  AL  =  0. 01,  A  =  1. 6E-03) 
PARAMETER  (ER  =  2.  2  ,  W  =  0. 711E-03) 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711,  1.0  OR  1.25  mm 

*  NOTE  :  PARAMETER  ER  CAN  BE  2.  2  OR  9.  8 

DO  15  FREQ  =  12.9075,12, 9125,3. 34E-4 
KK  =  CMPLX  (  0. 0,0. 0  ) 

*  START  SUMMATION  IN  ALPA  FROM  N=-30  TO  N=+30 
DO  25  L  =  1,61 

N  =  L  -  31 

ALPA  =  (N-0.5)*PI/A 

Kll  =  CMPLX  (  0.  0,0,  0  ) 

*  START  SUMMATION  IN  BETA  FROM  M=-3000  TO  M=+3000 

DO  35  I  =  1,6001 

M  =  I  -  3001 

BETA  =  (M-0.  5)*PI/(5.*AL) 

Kll  =  Kll  +  F1(ALPA,BETA,FREQ) 

35  CONTINUE 

KK  =  KK  +  Kll 
25  CONTINUE 

MAGKK  =  CABS((KK)) 

WRITE(6,100)  FREQ,  MAGKK 
100  FORMAT( '  ' ,3X,F25. 12,3X,F35.  10) 

15  CONTINUE 

STOP 
END 


*  FUNCTION  SUBPROGRAM  FI 

*  —  -  , 


COMPLEX  FUNCTION  Fl( ALPA, BETA, FREQ) 

COMPLEX  ZE , ZH , DE , NE , NH , COTl , COT2 , COT3 , GMl , GM2 , GM3 
COMPLEX  YE , YH , DHl , DH2 , DH , ZZZ 1 , ZZZ2 , ZZZ , JZ 1 

REAL  G 1 , G2 , G3 , AK 1 , AK2 , AKK 1 , AKK2 , GMM 1 , GMM2 , GMM3 , COTT 1 , COTTZ , COTTS 
REAL  YYE,YYH,ALPA,BETA,FREQ,AL,H,D,T,W,A,MO,PI,ER,EO 
REAL  X,Y,X1,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 
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PARAMETER  (PI  =  3.141592654,  EO  =  8. 85E-12) 

PARAMETER  (H  =  0. 66E-03,  T  =  0. 255E-03,  D  =  0. 66E-03) 
PARAMETER  (MO  =  PI*4.E-7,  AL  =  0. 01,  A  =  1.6E-03) 
PARAMETER  (W  =  0. 711E-03,  ER  =  2.  2  ) 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711,  1.0  OR  1.25  mm 

*  NOTE  :  PARAMETER  ER  CAN  BE  2.  2  OR  9. 8 

YYE  =  -1. /(2. E+9*PI*FREQ*E0) 

YE  =  CMPLX  (  0.0, YYE  ) 

YYH  =  2. E+9*PI*M0*FREQ 
YH  =  CMPLX(0.  0,YYH) 

AKl  =  M0*E0*((2.E+9*PI*FREQ)**2) 

AK2  =  M0*E0*ER*((2.E+9*PI*FREQ)**2) 

AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX(0.  0,G1) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX(0.  0,C0'm) 

C0T3  =  CMPLX(0. 0,C0TT3) 

ELSE 

GMl  =  CMPLX(G1,0.  0) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  1. /TANH(GMMl) 

C0TT3  =  1. /TANH(GMM3) 

COTl  =  CMPLX(COm,0.  0) 

C0T3  =  CMPLX(C0TT3,0.  0) 

END  IF 

G2  =  SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

C0TT2  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.  0,COTT2  ) 

ELSE 

GM2  =  CMPLX(G2,0.  0) 

GMM2  =  G2*T 

C0TT2  =  1. /TANH(GMM2) 

C0T2  =  CMPLX  (  C0TT2,0.  0  ) 

END  IF 

DE  =  C0T2*C0T3  +  C0T1*C0T3*(GM2/GM1)/ER  +  C0T1*C0T2  + 
NE  =  ((GM2/ER)*C0T3  +  GM3*COT2)*YE 
NH  =  GM2*C0T2  +  GM3*C0T3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 
DH2  =  GM2*GM3*COT2*COT3  +  GM2**2 
DH  =  DHl  +  DH2 
ZH  =  (  NH  /  DH  )*YH 
ZE  =  NE  /  DE 


(GM3/GM2)*ER 
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XI  =  -l./(ALPA**2  +  BETA**2) 

ZZZl  =  CMPLX  (  XI, 0.0  ) 

ZZZ2  =  (BETA**2)*ZE  +  (ALPA**2)*ZH 
ZZZ  =  ZZZ1*ZZZ2 
X  =  ALPA*W/2. 

Jll  =  2.*SIN(X)/X 
J12  =  3. /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  (2.*(1. -C0S(X)))/(X’'^2) 

J1  =  Jll  +  J12*(J13  -  J14  +  J15) 

Y  =  BETA*AL  /  2. 

J21  =  PI  *  COS(Y) 

J22  =  (PI/2. )**2  -  Y**2 

J2  =  J21  /  J22 

JJl  =  J1*J2 

JZl  =  CMPLX(JJ1,0.  0) 

FI  =  ZZZ*JZ1*JZ1 

RETURN 

END 
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B.  SECOND  ORDER  APPROXIMATION 


*  THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  THE  VALUE 

*  OF  DETERMINANT  OF  GREEN  FUNCTION  IN  SECOND  ORDER 

*  APPROXIMATION  WITH  W  =  0. 711  mm  and  ER  =  2.  2 

*  ■  ,  ■  , — ---  -  —  ■  -  .  -  -  --  -- 


*  MAIN  PROGRAM  : 

*  -  ■ — 


*  FKALPA, BETA, FREQ)  =  ZZZ  *  JZl  **  2 

*  F2(ALPA, BETA, FREQ)  =  ZZX  *  JZl  *  JXl 

*  F3(ALPA, BETA, FREQ)  =  ZXZ  *  JZl  *  JXl 

*  F4(ALPA, BETA, FREQ)  =  ZXX  *  JXl  **  2 

COMPLEX  K1 1 , K22 , K33 , K44 , KKl , KK2 , KK3 ,KK4 
COMPLEX  KKK1,KKK2,KK 
COMPLEX  F1,F2,F3,F4 
INTEGER  M,N,I,L 

REAL  MAGKK,A,W,D,T,H,MO,EO,ER,PI,AL,ALPA,BETA,FREQ 
PARAMETER  (  PI=3. 141592654,  E0=8.  85E-12  ) 

PARAMETER  (  H=0.66E-03,  T=0. 255E-03,  D=0. 66E-03  ) 
PARAMETER  (  M0=PI*4.  E-7,  AL=0.  01,  A  =1.  6E-03  ) 
PARAMETER  (  ER  =  2.  2  ,  W  =  0. 711E-03  ) 

*  NOTE  :  PAREMETER  W  CAN  BE  0.711,  1.0,  OR  1.25  mm 

*  NOTE  ;  PAREMETER  ER  CAN  BE  2.  2  OR  9.  8 

DO  15  FREQ  =  12. 8582191,12. 8582687,2. 48E-6 
KKl  =  CMPLXC  0. 0,0. 0  ) 

KK2  =  CMPLX(  0. 0,0. 0  ) 

KK3  =  CMPLXC  0. 0,0. 0  ) 

KK4  =  CMPLXC  0.  0,0.  0  ) 

*  START  SUMMATION  IN  ALPA  FROM  N=-30  TO  N=+30 
DO  25  L  =  1,  61 

N  =  L  -  31 

ALPA  =  CN-0. 5)*PI/A 

Kll  =  CMPLX  C  0. 0,0. 0  ) 

K22  =  CMPLX  C  0.0,0.  0  ) 

K33  =  CMPLX  C  0. 0,0. 0  ) 

K44  =  CMPLX  C  0. 0,0. 0  ) 

*  START  SUMMATION  IN  BETA  FROM  M=-3000  TO  M=+3000 

DO  35  I  =  1,  6001 

M  =  I  -  3001 

BETA  =  CM-0. 5)*PI/C5.*AL) 

Kll  =  Kll  +  FI C ALPA, BETA, FREQ) 

K22  =  K22  +  F2C ALPA, BETA, FREQ) 

K33  =  K33  +  F3C ALPA, BETA, FREQ) 

K44  =  K44  +  F4C ALPA, BETA, FREQ) 

35  CONTINUE 

KKl  =  KKl  +  Kll 

KK2  =  KK2  +  K22 

KK3  =  KK3  +  K33 

KK4  =  KK4  +  K44 

25  CONTINUE 
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KKKl  =  KK1*KK4 
KKK2  =  KK2*KK3 
KK  =  KKKl  -  KKK2 
MAGKK  =  CABS((KK)) 

WRITE(6,100)  FREQ,  MAGKK 
100  FORMATC^  ' ,2X,F16. 12,3X,F20. 10) 
15  CONTINUE 
STOP 
END 


FUNCTIONS  SUBPROGRAM  FI,  F2,  F3,  F4 


COMPLEX  FUNCTION  Fl( ALPA,BETA,FR£Q) 

COMPLEX  ZE , ZH , DE , NE , NH , COTl , COT2 , COT3 
COMPLEX  GM 1 , GM2 , GM3 , YE , YH , DH 1 , DH2 , DH 
COMPLEX  ZZZ1,ZZZ2,ZZZ,JZ1 

REAL  Gl,G2,G3,AKl,AK2,AKKl,AKK2,GMMl,GMM2,GMM3,COTTl,COTr2,COTT3 
REAL  YYE,YYH,ALPA,BETA,FREQ,AL,H,T,D,W,A,MO,PI,ER,EO 
REAL  X,Y,X1,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 
PARAMETER  (  PI=3.  141592654,  ER=2.  2,  E0=8.85E-12  ) 

PARAMETER  (  H=0.  66E-03,  T=0.255E-03,  D=0.  66E-03  ) 

PARAMETER  (  M0=PI*4.E-7,  AL=0. 01  ) 

PARAMETER  (  A=1.6E-03,  W=0.  711E-03  ) 

NOTE  :  PAREMETER  W  CAN  BE  0.711,  1.0,  OR  1.25  mm 
NOTE  :  PAREMETER  ER  CAN  BE  2.  2  OR  9.  8 

YYE  =  -1. /  (2.E+9*PI*FREQ*E0) 

YE  =  CMPLX  (0.0, YYE) 

YYH  =  2. E+9*PI*M0*FREQ 
YH  =  CMPLX  (0.0, YYH) 

AKl  =  M0*E0*((2.E+9*PI*FREQ)**2) 

AK2  =  M0*E0*ER*((2.  E+9*PI*FREQ)**2) 

AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0. 0,G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

COTTl  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX(0.  0,C0Tn) 

COT3  =  CMPLX(0.  0,C0TT3) 

ELSE 

GMl  =  CMPLX  (  G1,0. 0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

COTTl  =1./  TANH(GMMl) 

COTT3  =  1./  TANH(GMM3) 

COTl  =  CMPLX  (  COTTl, 0.0  ) 

COT3  =  CMPLX  (  COTT3,0.  0  ) 
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END  IF 

G2=SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

C0Tr2  =  -1. /TAN(GMM2) 

COT2  =  CMPLX  (  0.  0,COTT2  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

COTT2  =  1.  /  TANH(GMM2) 

C0T2  =  CMPLX  (  C0TT2,0.  0  ) 

END  IF 

DE  =  C0T2*C0T3+C0T1*C0T3*(GM2/GM1)/ER+C0T1*C0T2+(GM3/GM2)*ER 
NE  =  ((GM2/ER)*COT3+GM3*COT2)*YE 
NH=  GM2*COT2  +  GM3*COT3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 

DH2  =  GM2*GM3*COT2*COT3  +  GM2**2 

DH  =  DH1+DH2 

ZH  =  (NH  /  DH)  *  YH 

ZE  =  NE  /  DE 

XI  =  -1. /(  ALPA**2  +  BETA**2  ) 

ZZZl  =  CMPLX  (  XI, 0. 0  ) 

ZZZ2  =  (BETA**2)*ZE  +  (ALPA**2)*ZH 
ZZZ  =  ZZZ1*ZZZ2 
X  =  ALPA*W/2. 

Y  =  BETA*AL  /  2. 

Jll  =  2.*SIN(X)/(X) 

J12  =  3. /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  (2.*(1. -COS(X)))  /  (X**2) 

J1  =  J11+J12*(J13-J14+J15) 

J21  =  PI*COS(Y) 

J22  =  (PI/2. )**2  -  (Y**2) 

J2  =  J21/J22 
JJl  =  J1*J2 

JZl  =  CMPLX  (  JJl, 0.0  ) 

FI  =  ZZZ* JZl* JZl 

RETURN 

END 

COMPLEX  FUNCTION  F2(ALPA, BETA, FREQ) 

COMPLEX  ZE , ZH , DE , NE , NH , COTl , COT2 , COT3 
COMPLEX  GM 1 , GM2 , GM3 , YE , YH , DH 1 , DH2 , DH 
COMPLEX  ZZX,JZ1,JX1 

REAL  Gl,G2,G3,AKl,AK2,AKKl,AKK2,GMMl,GMM2,GMM3,COTTl,COTT2,COTr3 

REAL  YYE,YYH,ALPA,BETA,FREQ,AL,H,T,D,W,A,MO,PI,ER,EO 

REAL  X,Y,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 

REAL  J4 1 , J42 , J4 , J3 1 , J32 , J3 , JJ4 , ZZXl 

PARAMETER  (  PI=3. 141592654,  ER=2. 2,  E0=8. 85E-12  ) 

PARAMETER  (  H=0.66E-03,  T=0.  255E-03,  D=0.  66E-03  ) 

PARAMETER  (  M0=PI*4.E-7,  AL=0.  01  ) 

PARAMETER  (  A=0. 16E-02,  W=0.  711E-03  ) 

*  NOTE  :  PAREMETER  W  CAN  BE  0.711,  1.0,  OR  1.2)  mm 

*  NOTE  :  PAREMETER  ER  CAN  BE  2.  2  OR  9.  8 
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YYE  =  -1. /(2.E+9*PI*FREQ*E0) 

YE  =  CMPLX  (0.  O.YYE) 

YYH  =  2.  E+9*PI*M0*FREQ 
YH  =  CMPLX  (0. O.YYH) 

AKl  =  M0*E0*((2.  E+9*PI*FREQ)**2) 

AK2  =  M0*E0*ER*((2.  E+9*PI*FREQ)**2) 

AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0. 0,G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX(0.  O.COTTl) 

C0T3  =  CMPLX(0.  0,C0TT3) 

ELSE 

GMl  =  CMPLX  (  G1,0. 0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  1. /TANH(GMMl) 

C0TT3  =  1. /TANH(GMM3) 

COTl  =  CMPLX  (  com, 0.0  ) 

C0T3  =  CMPLX  (  COTT3,0.0  ) 

END  IF 

G2=SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

C0TT2  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.  0,C0TT2  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

C0Tr2  =  1. /TANH(GMM2) 

C0T2  =  CMPLX  (  C0TT2,0.  0  ) 

END  IF 

DE  =  C0T2*C0T3+C0T1*C0T3*(GM2/GM1)/ER+C0T1*C0T2+(GM3/GM2)*ER 
NE  =  ((GM2/ER)*COT3+GM3*COT2)*YE 
NH  =  GM2*C0T2  +  GM3*C0T3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 
DH2  =  GM2*GM3*COT2*COT3  +  GM2**2 
DH  =  DH1+DH2 
ZH  =  (NH/DH)*YH 
ZE  =  NE/DE 

ZZXl  =  -l.*(ALPA*BETA)/(ALPA**2+BETA**2) 

ZZX  =  ZZX1*(ZE-ZH) 

X  =  ALPA*W/2. 

Y  =  BETA*AL/2. 

J31  =  2.*PI*SIN(X) 

J32  =  X**2  -  PI**2 
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J3  =  J31/J32 
J41  =  COS(Y)/(Y) 

J42  =  SIN(Y)/(Y**2) 

J4  =  J41-J42 
JJ4  =  J3*J4 

JXl  =  CMPLX  (  JJ4,0.  0  ) 

Jll  =  2.*SIN(X)/(X) 

J12  =  3.  /(X**2) 

J13  =  COS(X) 

J14  =  2.*SIN(X)/(X) 

J15  =  (2.*(1. -C0S(X)))/(X**2) 

J1  =  J11+J12*(J13-J14+J15) 

J21  =  PI*COS(Y) 

J22  =  (PI/2)**2-Y**2 

J2  =  J21/J22 

JJl  =  J1*J2 

JZl  =  CMPLX  (JJl.O.  0  ) 

F2  =  ZZX*JX1*JZ1 

RETURN 

END 

COMPLEX  FUNCTION  F3(ALPA, BETA, FREQ) 

COMPLEX  ZE,ZH,DE,NE,NH,C0T1,C0T2,C0T3 
COMPLEX  GM 1 , GM2 , GM3 , YE , YH , DH 1 , DH2 , DH 
COMPLEX  ZXZ,JZ1,JX1 

REAL  Gl,G2,G3,AKl,AK2,AKKl,AKK2,GMMl,GMM2,GMM3,COTTl,COTT2,COTT3 

REAL  YYE,YYH,ALPA,BETA,FREQ,AL,H,T,D,W,A,MO,PI,ER,EO 

REAL  X,Y,J11,J12,J13,J14,J15,J1,J21,J22,J2,JJ1 

REAL  J41,J42,J4,J31,J32,J3,JJ4,ZXZ1 

PARAMETER  (  PI=3. 141592654,  ER=2.  2,  E0=8.85E-12  ) 

PARAMETER  (  H=0.  66E-03.  T=0.  255E-03,  D=0. 66E-03  ) 

PARAMETER  (  M0=PI*4.E-7,  AL=0.01  ) 

PARAMETER  (  A=1.6E-03,  W=0.  711E-03  ) 

*  NOTE  :  PAREMETER  W  CAN  BE  0.711,  1.0,  OR  1.25  mm 

*  NOTE  :  PAREMETER  ER  CAN  BE  2.  2  OR  9.  8 

YYE  =  -l./(2.E+9*PI*FREQ*E0) 

YE  =  CMPLX  (0.  0,YYE) 

YYH  =  2. E+9*PI*M0*FREQ 
YH  =  CMPLX  (0. 0,YYH) 

AKl  =  M0*E0*((2.  E+9*PI*FREQ)**2) 

AK2  =  M0*E0*ER*((2.E+9*PI*FREQ)**2) 

AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA’^2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0.0, G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX(0.  0,C0Tn) 

C0T3  =  CMPLXCO.  0,00713) 

ELSE 
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GMl  =  CMPLX  (  G1,0.  0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  1. /TANH(GMMl) 

C0TT3  =  1.  /TANH(GMM3) 
coTi  =  CMPLX  (  com, 0.0  ) 

C0T3  =  CMPLX  (  C0TT3,0.  0  ) 

END  IF 

G2=SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

com  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.0, com  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

com  =  1.  /TANH(GMM2) 

C0T2  =  CMPLX  (  com, 0.0  ) 

END  IF 

DE  =  COT2*COT3+COT1*COT3*(GM2/GM1)/ER+COT1*COT2+(GM3/GM2)*ER 
NE  =  ((GM2/ER)*C0T3+GM3*C0T2)*YE 
NH  =  GM2*C0T2  +  GM3*COT3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 

DH2  =  GM2*GM3*C0T2*C0T3  +  GM2’'«^2 

DH  =  DH1+DH2 

ZH  =  (NH  /  DH)  *  YH 

ZL  =  NE  /  DE 

ZXZl  =  -l.*(ALPA*BETA)/(ALPA**2  +  BETA**2) 

ZXZ  =  ZXZ1*(ZE-ZH) 

X  =  ALPA*W/2. 
y  =  BETA*AL  /  2. 

Jll  =  2. *SIN(X)/(X) 

J12  =  3. /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  (2. 0*(1. O-COS(X)))  /  (X**2) 

J1  =  J11+J12*(J13-J14+J15) 

J21  =  PI*COS(Y) 

J22  =  (PI/2. )**2  -  (Y**2) 

J2  =  J21/J22 
JJl  =  J1*J2 

JZl  =  CMPLX  (  JJl, 0.0  ) 

J31  =  2.*PI*SIN(X) 

J32  =  X**2  -  PI**2 
J3  =  J31  /  J32 
J41  =  COS(Y)/(Y) 

J42  =  SIN(Y)/(Y**2) 

J4  =  J41-J42 
JJ4  =  J3*J4 

JXl  =  CMPLX  (  JJ4,0.0  ) 

F3  =  ZXZ*JZ1*JX1 

RETURN 

END 
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COMPLEX  FUNCTION  F4(ALPA, BETA, FREQ) 

COMPLEX  ZE , ZH , DE , NE , NH , COT 1 , COT2 , COT3 
COMPLEX  GMl , GM2 , GM3 , YE , YH , DH 1 , DH2 , DH 
COMPLEX  ZXX1,ZXX2,ZXX,JX1 

REAL  G 1 , G2 , G3 , AKl , AK2 , AKK 1 , AKK2 , GMMl , GMM2 , GMM3 , COTTl , COTT2 , COTTS 
REAL  YYE,YYH,ALPA,BETA,FREQ,AL,H,T,D,W,A,MO,PI,ER,EO 
REAL  X,Y,X1,J31,J32,J3,J41,J42,J4,JJ4 
PARAMETER  (  PI=3. 141592654,  ER=2.  2,  E0=8.  85E-12  ) 

PARAMETER  (  H=0.  66E-03,  T=0. 255E-03,  D=0.  66E-03  ) 

PARAMETER  (  M0=PI*4.  E-7,  AL=0.  01  ) 

PARAMETER  (  A=l.  6E-03,  W=0. 711E-03  ) 

*  NOTE  :  PAREMETER  W  CAN  BE  0.711,  1.0,  OR  1.25  mm 

*  NOTE  :  PAREMETER  ER  CAN  BE  2.  2  OR  9.  8 

YYE  =  -1. /(2.E+9*PI*FREQ*E0) 

YE  =  CMPLX  (0.  0,YYE) 

YYH  =  2. E+9*PI*M0*FREQ 
YH  =  CMPLX  (0.  0,YYH) 

AKl  =  M0*E0*((2.E+9*PI*FREQ)**2) 

AK2  =  M0*E0*ER*((2.  E+9*PI*FREQ)**2) 

AKKl  =  ALPA**2  +  !i£TA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0.0, G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

COTTl  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLXCO. 0, COTTl) 

COT3  =  CMPLXCO.  0,C0TT3) 

ELSE 

GMl  =  CMPLX  (  G1,0. 0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

COTTl  =  1. /TANH(GMMl) 

COTr3  =  1. /TANH(GMM3) 

COTl  =  CMPLX  (  COTTl, 0.0  ) 

COT3  =  CMPLX  (  cons, 0.0  ) 

ENDIF 

G2=SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

COTr2  =  -1. /TAN(GMM2) 

COT2  =  CMPLX  (  0. 0,C0TT2  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

COTT2  =  1. /TANH(GMM2) 

COT2  =  CMPLX  (  C0TT2,0.  0  ) 

END  IF 

DE  =  COT2*COT3+COT1*COT3*(GM2/GM1)/ER+COT1*COT2+(GM3/GM2)*ER 
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NE  =  ((GM2/ER)*C0T3+GM3*C0T2)*yE 
NH  =  GM2*C0T2  +  GM3*C0T3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 

DH2  =  GM2*GM3*C0T2*C0T3  +  GM2**2 

DH  =  DH1+DH2 

ZH  =  (NH  /  DH)  *  YH 

7F  =  NF  /  DF 

Xi  =  -1. /(ALPA**2  +  BETA**2) 

ZXXl  =  CMPLX  (XI, 0.0) 

ZXX2  =  (ALPA**2)*ZE  +  (BETA**2)*ZH 
ZXX  =  ZXX1*ZXX2 
X  =  ALPA*W/2. 

Y  =  BETA*AL/2. 

J31  =  2.*PI*SIN(X) 

J32  =  X**2  -  PI**2 
J3  =  J31  /  J32 
J41  =  COS(Y)/(Y) 

J42  =  SIN(Y)/(Y**2) 

J4  =  J41-J42 
JJ4  =  J3*J4 

JXl  =  CMPLX  (  JJ4,0.  0  ) 

F4  =  ZXX* JXl* JXl 

RETURN 

END 
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APPENDIX  C.  PROGRAMS  FOR  CALCULATION  OF 
PROPAGATION  CONSTANT 


* 

* 

ie 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 

* 


VARIABLE  DEFINITION: 

I  =  REPRESENT  THE  REGION  1,2,3  IN  RESONATOR 

ALPA  =  DISCRETE  TRANSFORM  VARIABLE  FOR  DOMINANT  MODE  (x-direction) 

BETA  =  DISCRETE  TRANSFORM  VARIABLE  FOR  DOMINANT  MODE  ( z-direction) 

A  =  A  HALF  WIDTH  OF  RESONATOR  DEVICE  =  1.6  tnm 

D  =  HEIGHT  OF  RESONATOR  BELOW  SUBSTRATE  =  0.  66  mm 

T  =  THE  THICKNESS  OF  SUBSTRATE  =  0. 255  mm 

W  =  THE  WIDTH  OF  STRIP,  CAN  BE  0.711  mm,  1.0  mm,  or  1.25  mm 

H  =  THE  HEIGHT  OF  RESONATOR  ABOVE  SUBSTRATE  =  0.66  mm 

AL  =  THE  LENGTH  OF  STRIPLINE  ,  ASSUME  =  10  mm 

AKI  =  RADIAN  IREQUENCY  *  SQRT  OF  (PERMimVITY*PERMEABILITY) 

MO  =  PERMEABILITY  IN  VACCUM  =  PI*4.E-7 

MR  =  PERMEABILITY  IN  A  SUBSTRATE  ASSUME  =  1 

EO  =  PERMITTIVITY  IN  VACCUM  =  (l/(36*PI))*E-9  =  8. 85E-12 

ER  =  PERMITIVITY  IN  A  SUBSTRATE,  CAN  BE  2.2  OR  9.8 

FREQ  =  FREQUENCY  OPERATION  IN  GHZ 

DE  =  DENOMINATOR  OF  ZE 

DE  =  NUMERATOR  OF  ZE 

DH1,DH2,DH  =  DENOMINATOR  OF  ZH 

NH  =  NUMERATOR  OF  ZH 

ZE  =  IMPEDANCE  OF  ELECTRIC  FIELD,  eqn  (Ig) 

ZH  =  IMPEDANCE  OF  MAGNETIC  FIELD,  eqn  (Ih) 

ZZZ,ZZX,2XZ,AND  ZXX  ARE  THE  GREEN’S  IMPEDANCE  FUNCTION  (lc,d,e,f) 
COTl  =  CTl,  EQN  (li) 

C0T2  =  CT2,  EQN  (Ij) 

COT3  =  CT3,  EQN  (Ik) 

J1  =  EQN  (8a) 

J2  =  EQN  (8b) 

JJl  =  Jzl  EQN  (6a) 

JJ4  =  Jxl  EQN  (6b) 

JZl  =  BASIS  FUNCTION  OF  CURRENT  COMPONENT  IN  Z  -  DIRECTION 

JXl  =  BASIS  FUNCTION  OF  CURRENT  COMPONENT  IN  X  -  DIRECTION 

K11,K12,K21,K22  =  MATRIX  ELEMENT 

Jl,  J3  =  ASSUME  CURRENT  DISTRIBUTION  IN  X  -  DIRECTION 
J2,  J4  =  ASSUME  CURRENT  DISTRIBUTION  IN  Z  -  DIRECTION 
ALL  COMPUTATION  IN  STANDARD  "MKS" 
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A.  FIRST  ORDER  APPROXIMATION 


*  THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  CONSTANT  PROPAGATION 

*  IN  FIRST  ORDER  APPROXIMATION  WITH  W  =  0. 711  AND  ER  =  2.  2 


*  MAIN  PROGRAM  : 


*  FUNCTION  :  F1(ALPA, BETA, FREQ)  =  ZZZ  *  JZl  **  2 
COMPLEX  F1,K11 

REAL  B0,B00,MAGKK,ALPA,BETA,M0 
REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO 
INTEGER  N,L 

PARAMETER  (  PI=3.  141592654,  ER=2. 2,  E0=8. 85E-12,  H=0. 66E-03  ) 
PARAMETER  (  T=0.  255E-03,  D=0.  66E-03,  M0=PI*4.E-7,  A=0.  16E-02  ) 
PARAMFTFR  (  W=n  711F-3  \ 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711  MM,  1.0  MM,  1.25  MM 

*  NOTE  :  PARAMETER  ER  CAN  BE  2.2  OR  9.8 

FREQ  =  12.9095783 

OMG  =  FREQ*PI*2,  E+9 

BO  =  (0MG/3.E08)*SQRT((ER+l)/2) 

*  START  SWEEP  THE  BETA  VALUE  FROM  84%  TO  90% 

DO  20  BOO  =  0. 84,0. 9,1. 2E-3 

BETA  =  B0*B00 
Kll  =  ,0. 0,0.0) 

*  START  SUMMATION  IN  ALPA  FROM  N=-40  TO  N=+40 

DO  10  L  =  1,  81 

N  =  L  -  41 

ALPA  =  (N  -  0.5)*PI/A 

Kll  =  Kll  +  Fl( ALPA, BETA, FREQ) 

10  CONTINUE 

MAGKK  =  CABS((K11)) 

WRITE(6,100)  FREQ, BETA, MAGKK 
100  FORMATC  ' ,3X,F14. 9,3X,F18. 9,3X,F15. 9) 

20  CONTINUE 

STOP 
END 

*  FUNCTION  SUBPROGRAM  FI 


COMPLEX  FUNCTION  F1(ALPA,BETA,FREQ) 

COMPLEX  GMl , GM2 , GM3 , COTl , COT2 , COT3 , DE , ZE , DH , ZH , NE , NH , DHl , DH2 , YE 
COMPLEX  JZ1,ZZZ1,ZZZ2,ZZZ,YH 

REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO,Gl,G2,G3,COm,COTT2,COTT3,YYE 

REAL  YYH , GMM 1 , GMM2 , GMM3 , ALPA , AKK 1 , AKK2 , X 1 , X 

REAL  J11,J12,J13,J14,J15,J1,M0,AK1,AK2,BETA 

PARAMETER  (  PI=3. 141592654,  ER=2. 2,  E0=8.85E-12,  H=0. 66E-03  ) 

PARAMETER  (  T=0.255E-03,  D=0.66E-03,  M0=PI*4.E-7,  A=0. 16E-02  ) 
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★ 

* 


PARAMFTFR  (  U=n  711F-T 

NOTE  :  PARAMETER  W  CAN  BE  0.711  MM,  1.0  MM,  1.25  MM 
NOTE  :  PARAMETER  ER  CAN  BE  2.  2  OR  9.  8 


<1* 


OMG  =  FREQ*PI*2.  E+9 
YYE  =  -1. /(OMG*EO) 

YE  =  CMPLX  (0.0, YYE) 

YYH  =  OMG*MO 

YH  =  CMPLX  (0.  0,YYH) 

AKl  =  (OMG**2)*MO*EO 
AK2  =  (OMG**2)*MO*EO*ER 
AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0.0, G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX  (  0.  0,C0TT1  ) 

C0T3  =  CMPLX  (  0.0,COTT3  ) 

ELSE 

GMl  =  CMPLX  (  G1,0.  0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G1*D 

com  =  1. /TANH(GMMl) 

com  =  1. /TANH(GMM3) 

COTl  =  CMPLX  (  com, 0.0  ) 

C0T3  =  CMPLX  (  com, 0.0  ) 

END  IF 

G2  =  SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

com  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.0, com  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

com  =  1. /TANH(GMM2) 

C0T2  =  CMPLX  (  com, 0.0  ) 

END  IF 

DE  =  C0T2*C0T3+C0T1*C0T3*(GM2/GM1)/ER+C0T1*C0T2+(GM3/GM2)*ER 
NE  =  (GM2*C0T3/ER+GM3*C0T2)*YE 
ZE  =  NE/DE 

NH  =  GM2*C0T2  +  GM3*C0T3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 

DH2  =  GM2*GM3*C0T2*C0T3  +  GM2**2 

DH  =  DHl  +  DH2 

ZH  =  (NH/DH)*YH 

XI  =  -1.  /  (ALPA**2  +  BETA**2) 

ZZZl  =  CMPLX  (XI, 0.  0) 
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ZZZ2  =  (BETA**2)*ZE  +  (ALPA**2)*ZH 
ZZZ  =  ZZZ1*ZZZ2 
X  =  ALPA*W/2. 

Jll  =  2.*SIN(X)/(X) 

J12  =  3. /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  2.*(1.0-C0S(X))/(X**2) 

J1  =  J11+J12*(J13-J14+J15) 

JZl  =  CMPLX  (  J1,0. 0  ) 

FI  =  ZZZ*JZ1**2 

RETURN 

END 
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B.  SECOND  ORDER  APPROXIMATION 


*  THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  CONSTANT  PROPAGATION 

*  IN  SECOND  ORDER  APPROXIMATION  WITH  W  =  0. 711  AND  ER  =  2.  2 

*  i  ■■  - '  -  '■  —  —  --  - — - 


*  MAIN  PROGRAM  : 

*  . . -  - 


* 

* 

* 

* 


F1(ALPA,BETA,FREQ) 
F2(ALPA, BETA, FREQ) 
F3(ALPA, BETA, FREQ) 
F4( ALP A, BETA, FREQ) 


=  zzz 

* 

JZl 

**  2 

=  zzx 

* 

JZl 

*  JXl 

=  zxz 

* 

JZl 

*  JXl 

=  zxx 

* 

JXl 

**  2 

COMPLEX  F1,F2,F3,F4 
COMPLEX  KK,K11,K22,K33,K44 
COMPLEX  KKK1,KKK2 
REAL  MAGKK,BOO,BO,ALPA,BETA,MO 
REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO 
INTEGER  N,L 

PARAMETER  (  PI=3. 141592654,  E0=8.  85E-12,  H=0.  66E-03  ) 

PARAMETER  (  T=0.255E-03,  D=0.  66E-03,  M0=PI*4.  E-7,  A=0.  16E-02  ) 
PARAMETER  (  W=0. 711  ,ER=2. 2  ) 

NOTE  :  PARAMETER  W  CAN  BE  0.711  ,  1. 0  ,  OR  1.25 
NOTE  :  PARAMETER  ER  CAN  BE  2.  2  OR  9.8 


FREQ  =  12. 858241 

OMG  =  FREQ*PI*2.  E+9 

BO  =  COMG/3.E08)*SQRT((ER+l)/2) 

*  START  SWEEP  THE  BETA  VALUE  FROM  88.2%  TO  89.7% 
DO  20  BOO  =  0. 882,0. 897,7. 5E-6 

BETA  =  B0*B00 
Kll  =  (0.  0,0.  0) 

K22  =  (0. 0,0. 0) 

K33  =  (0.0,0.  0) 

K44  =  (0. 0,0. 0) 

*  START  SUMMATION  IN  ALPA  FROM  N=-40  TO  N=+40 

DO  10  L  =  1,  81 

N  =  L  -  41 

ALPA  =  (N  -  0. 5)*PI/A 

Kll  =  Kll  +  FKALPA, BETA, FREQ) 

K22  =  K22  +  F2( ALPA, BETA, FREQ) 

K33  =  K33  +  F3( ALPA, BETA, FREQ) 

K44  =  K44  +  F4( ALPA, BETA, FREQ) 

10  CONTINUE 

KKKl  =  K11*K44 
KKK2  =  K22*K33 
KK  =  KKK1-KKK2 
MAGKK  =  CABS((KK)) 

WRITE(6,100)  FREQ, BETA, MAGKK 
100  FORMATC  ' ,3X,F14. 9,3X,F18.  9,3X,F15. 9) 

20  CONTINUE 
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STOP 

END 


*  FUNCTIONS  SUBPROGRAM  FI,  F2,  F3,  F4 

*  - - 


COMPLEX  FUNCTION  F1(ALPA, BETA, FREQ) 

COMPLEX  GMl , GM2 , GM3 , COTl , COT2 , COT3 , DE , ZE , DH , ZH , NE , NH , DH 1 , DH2 , YE 
COMPLEX  YH,J21,ZZZ1,ZZZ2,ZZZ 

REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO,Gl,G2,G3,COm,COTr2,COTT3,YYE 
REAL  YYH,GMM1,GMM2,GMM3,ALPA,AKK1,AKK2,X1,X 
REAL  JH,J12,J13,J14,J15,J1,M0,AK1,AK2,BETA 
PARAMETER  (  PI=3.  141592654,  E0=8.85E-12,  H=0.66E-03  ) 

PARAMETER  (  T=0.255E-03,  r)=0.66E-03,  M0=PI*4.E-7,  A=0.  16E-02  ) 
PARAMETER  (  W=0.  711  ,ER=2. 2  ) 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711  ,  1. 0  ,  OR  1.25 

*  NOTE  :  PARAMETER  ER  CAN  BE  2. 2  OR  9.  8 


OMG  =  FREQ*PI*2.  E+9 
YYE  =  -l./(OMG*EO) 

YE  =  CMPLX  (0.0, YYE) 

YYH  =  OMG*MO 

YH  =  CMPLX  (0. 0,YYH) 

AKl  =  (0MG**2)*M0*E0 
AK2  =  (OMG**2)*MO*EO*ER 
AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0.0, G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  -1.  /TAN(GMMl) 

COTT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX  (  0.0, com  ) 

C0T3  =  CMPLX  (  0. 0,C0TT3  ) 

ELSE 

GMl  =  CMPLX  (  G1,0. 0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G1*D 

com  =  1. /TANH(GMMl) 

C0TT3  =  1. /TANH(GMM3) 

COTl  =  CMPLX  (  com, 0.0  ) 

C0T3  =  CMPLX  (  com, 0.0  ) 

END  IF 

G2  =  SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

C0TT2  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.0, com  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 
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GMM2  =  G2*T 

com  =  1. /TANH(GMM2) 

COT2  =  CMPLX  (  com, 0.0  ) 

END  IF 

DE  =  C0T2*C0T3+C0T1*C0T3*(GM2/GM1)/ER+C0T1*C0T2+(GM3/GM2)*ER 
NE  =  (GM2*C0T3/ER+GM3*C0T2)*YE 
ZE  =  NE/DE 

NH  =  GM2*C0T2  +  GM3*C0T3 

DHl  =  GM1*GM2*CCT1*C0T2  +  GM1*GM3*C0T1*C0T3 
DH2  =  GM2*GM3*C0T2*C0T3  +  GM2**2 
DH  =  DHl  +  DH2 
ZH  =  (NH/DH)*YH 

XI  =  -1.  /  (ALPA**2  +  BETA**2) 

ZZZl  =  CMPLX  (XI, 0.0) 

ZZZ2  =  (BETA**2)*ZE  +  (ALPA**2)*ZH 
ZZZ  =  ZZZ1*ZZZ2 
X  =  ALPA*W/2. 

Jll  =  2.*SIN(X)/(X) 

J12  =  3.  /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  2.*(1.  0-COS(X))/(X**2) 

J1  =  J11+J12*(J13-J14+J15) 

JZl  =  CMPLX  (  J1,0. 0  ) 

FI  =  ZZZ*JZ1**2 

RETURN 

END 

COMPLEX  FUNCTION  F2(ALPA, BETA, FREQ) 

COMPLEX  GM 1 , GM2 , GM3 , COTl , C0T2 , C0T3 , DE , ZE , DH , ZH , NE , NH , DHl , DH2 , YE 
COMPLEX  yH,JZl,JXl,ZZX 

REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO,Gl,G2,G3,COTn,COm,COTT3,YYE 
REAL  YYH , GMMl , GMM2 , GMM3 , ALPA , AKKl , AKK2 , X , ZZX 1 , BETA 
REAL  J31,J32,J3,J11,J12,J13,J14,J15,J1,M0,AK1,AK2 
PARAMETER  (  PI=3. 141592654,  E0=8.85E-12,  H=0. 66E-03  ) 

PARAMETER  (  T=0. 255E-03,  D=0.  66E-03,  M0=PI*4.E-7,  A=0.  16E-02  ) 
PARAMETER  (  W=0.  711  ,ER=2.  2  ) 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711  ,  1.  0  ,  OR  1.25 

*  NOTE  ;  PARAMETER  ER  CAN  BE  2. 2  OR  9. 8 

OMG  =  FREQ*PI*2.  E+9 
YYE  =  -l./(OMG*EO) 

YE  =  CMPLX  (0.  0,YYE) 

YYH  =  OMG*MO 

YH  =  CMPLX  (0.  0,YYH) 

AKl  =  (0MG**2)*M0*E0 
AK2  =  (OMG**2)*MO*EO*ER 
AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0. 0,G1  ) 

GM3  =  GMl 
GMMl  =  G1*H 
GMM3  =  G3*D 
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COTTl  =  -1, /TAN(GMMl) 

COTT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX  (  0.0, COTTl  ) 

C0T3  =  CMPLX  (  0. 0,C0TT3  ) 

ELSE 

GMl  =  CMPLX  (  G1,0. 0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G1*D 

COTTl  =  l./TANH(GMMl) 

C0TT3  =  1. /TANH(GMM3) 

COTl  =  CMPLX  (  COTTl, 0.0  ) 

C0T3  =  CMPLX  (  COTT3,0.  0  ) 

END  IF 

G2  =  SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0. 0,G2  ) 

GMM2  =  G2*T 

COTT2  =  -1. /TAN(GMM2) 

COT2  =  CMPLX  (  0. 0,C0TT2  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

C0TT2  =  1. /TANH(GMM2) 

C0T2  =  CMPLX  (  COTT2,0.0  ) 

END  IF 

DE  =  COT2*COT3+COTl*COT3*(GM2/GMl)/ER+COTl*COT2+(GM3/GM2)*ER 
NE  =  (GM2*C0T3/ER+GM3*C0T2)*YE 
ZE  *  NE/DE 

NH  =  GM2*COT2  +  GM3*C0T3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 

DH2  =  GM2*GM3*COT2*COT3  +  GM2**2 

DH  =  DHl  +  DH2 

ZH  =  (NH/DH)*YH 

ZZXl  =  -l.*(ALPA*BETA)/(ALPA**2+BETA**2) 

ZZX  =  ZZX1*(ZE-ZH) 

X  =  ALPA*W/2. 

Jll  =  2.*SIN(X)/(X) 

J12  =  3. /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  2.*(1.0-COS(X))/(X**2) 

J1  =  J11+J12*(J13-J14+J15) 

J31  =  2.*PI*SIN(X) 

J32  =  X**2  -  PI**2 

J3  =  J31/J32 

JZl  =  CMPLX  (  J1,0.  0  ) 

JXl  =  CMPLX(J3,0.  0) 

F2  =  ZZX* JZl* JXl 

RETURN 

END 


COMPLEX  FUNCTION  F3( ALFA, BETA, FREQ) 

COMPLEX  GM 1 , GM2 , GM  3 , COT 1 , C0T2 , C0T3 , DE , ZE , DH , ZH , NE , NH , DH 1 , DH2 . YE 
COMPLEX  YH,JZ1,JX1,ZXZ 

REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO,Gl,G2,G3,COTTl,COTT2,COTT3,YyE 
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REAL  YYH , GMMl , GMM2 , GMM3 , ALFA , AKKl , AKK2 , X , BETA 
REAL  J11,J12,J13,J14,J15,J1,J31,J32,J3,M0,AK1,AK2 
PARAMETER  (  PI=3.  141592654,  E0=8. 85E-12,  H=0. 66E-03  ) 

PARAMETER  (  T=0. 255E-03,  D=0. 66E-03,  M0=PI*4.E-7,  A=0.  16E-02  ) 
PARAMETER  (  W=0.  711  ,ER=2.  2  ) 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711  ,  1. 0  ,  OR  1.25 

*  NOTE  :  PARAMETER  ER  CAN  BE  2.  2  OR  9.  8 

OMG  =  FREQ*PI*2.  E+9 
YYE  =  -l./(OMG*EO) 

YE  =  CMPLX  (0.0, YYE) 

YYH  =  OMG*MO 

YH  =  CMPLX  (0. 0,YYH) 

AKl  =  (OMG**2)*MO*EO 
AK2  =  (0MG**2)*M0*E0*ER 
AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABSCAKKl)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0.0, G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

COTTl  =  -1. /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 

COTl  =  CMPLX  (  0.0, COTTl  ) 

C0T3  =  CMPLX  (  0. o,corr3  ) 

ELSE 

GMl  =  CMPLX  (  G1,0.0  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G1*D 

COTTl  =  1. /TANH(GMMl) 

C0TT3  =  1.  /TANH(GMM3) 

COTl  =  CMPLX  (  COTTl, 0.0  ) 

C0T3  =  CMPLX  (  COTT3,0.  0  ) 

END  IF 

G2  =  SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

C0TT2  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.0,CO'IT2  ) 

ELSE 

GM2  =  CMPLX  (  G2,0. 0  ) 

GMM2  =  G2*T 

C0Tr2  =  1. /TANH(GMM2) 

C0T2  =  CMPLX  (  C0Tr2,0.  0  ) 

END  IF 

DE  =  C0T2*C0T3+C0Tl*C0T3*(GM2/GMl)/ER+r0Tl*C0T2+(GM3/GM2)*ER 
NE  =  (GM2*C0T3/ER+GM3*C0T2)*YE 
ZE  =  NE/DE 

NH  =  GM2*C0T2  +  GM3*C0T3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 
DH2  =  GM2*GM3*C0T2*C0T3  +  GM2**2 
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<1* 


DH  =  DHl  +  DH2 
ZH  =  (NH/DH)*YH 

ZXZ  =  -l.*(ALPA*BETA)/(ALPA**2+BETA**2) 

ZXZ  =  ZXZ*(ZE-ZH) 

X  =  ALPA*W/2. 

J31  =  2.*PI*SIN(X) 

J32  =  X**2  -  PI**2 
J3  =  J31/J32 
JXl  =  CMPLX(J3,0.  0) 

Jll  =  2.*SIN(X)/(X) 

J12  =  3. /(X**2) 

J13  =  COS(X) 

J14  =  Jll 

J15  =  2.  *(1.  0-C0S(X))/(X**2) 

J1  =  J11+J12*(J13-J14+J15) 

JZl  =  CMPLX  (  J1,0.  0  ) 

F3  =  2XZ*JZ1*JX1 
RETURN 
END 

COMPLEX  FUNCTION  F4(ALPA, BETA, FREQ) 

COMPLEX  GM1,GM2,GM3,COT1,COT2,COT3,DE,ZZ,DH,ZH,NE,NH,DH1,DH2,YE 
COMPLEX  YH,JX1,ZXX1,ZXX2,ZXX 

REAL  A,D,T,H,PI,ER,W,FREQ,OMG,EO,Gl,G2,G3,COTTl,COrr2,COTT3,YYE 
REAL  YYH , GMM 1 , GMM2 , GMM3 , ALPA , AKK 1 , AKK2 , AK 1 , AK2 
REAL  J31,J32,J3,M0,X,X1,BETA 

PARAMETER  (  PI=3. 141592654 ,  E0=8. 85E-12,  H=0.  66E-03  ) 

PARAMETER  (  T=0.255E-03,  D=0.66E-03,  M0=PI*4.E-7,  A=0.  16E-02  ) 
PARAMETER  (  W=0. 711  ,ER=2. 2  ) 

*  NOTE  :  PARAMETER  W  CAN  BE  0.711  ,  1.  0  ,  OR  1.25 

*  NOTE  :  PARAMETER  ER  CAN  BE  2.  2  OR  9.  8 

OMG  =  FREQ*PI*2.  E+9 
YYE  =  -1. /(OMG*EO) 

YE  =  CMPLX  (0.0, YYE) 

YYH  =  OMG*MO 

YH  =  CMPLX  (0.  0,YYH) 

AKl  =  (OMG**2)*MO*EO 
AK2  =  (0MG**2)*M0*E0*ER 
AKKl  =  ALPA**2  +  BETA**2  -  AKl 
AKK2  =  ALPA**2  +  BETA**2  -  AK2 
G1  =  SQRT(ABS(AKK1)) 

G3  =  G1 

IF  (AKKl  .LT.  0.0)  THEN 
GMl  =  CMPLX  (  0. 0,G1  ) 

GM3  =  GMl 

GMMl  =  G1*H 

GMM3  =  G3*D 

com  =  -1.  /TAN(GMMl) 

C0TT3  =  -1. /TAN(GMM3) 
coTi  =  CMPLX  (  0.0, com  ) 

C0T3  =  CMPLX  (  0.  0,C0TT3  ) 

ELSE 

GMl  =  CMPLX  (  G1,0. 0  ) 

GM3  =  GMl 
GMMl  =  G1*H 
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GMM3  =  G1*D 

com  =  1.  /TANH(GMMl) 

COTT3  =  1. /TANH(GMM3) 
coTi  =  CMPLx  (  com, 0.0  ) 

C0T3  =  CMPLX  (  C0TT3,0.  0  ) 

END  IF 

G2  =  SQRT(ABS(AKK2)) 

IF  (AKK2  .LT.  0.0)  THEN 
GM2  =  CMPLX  (  0.0, G2  ) 

GMM2  =  G2*T 

C0TT2  =  -1. /TAN(GMM2) 

C0T2  =  CMPLX  (  0.  0,COTT2  ) 

ELSE 

GM2  =  CMPLX  (  G2,0.0  ) 

GMM2  =  G2*T 

com  =  1. /TANH(GMM2) 

C0T2  =  CMPLX  (  com, 0.0  ) 

END  IF 

DE  =  C0T2*C0T3+C0T1*C0T3*(GM2/GM1)/ER+C0T1*C0T2+(GM3/GM2)*ER 

NE  =  (GM2*C0T3/ER+GM3*C0T2)*YE 
ZE  =  NE/DE 

NH  =  GM2*C0T2  +  GM3*COT3 

DHl  =  GM1*GM2*C0T1*C0T2  +  GM1*GM3*C0T1*C0T3 

DH2  =  GM2*GM3*C0T2*C0T3  +  GM2**2 

DH  =  DHl  +  DH2 

ZH  =  (NH/DH)*YH 

XI  =  -l./(ALPA**2  +  BETA**2) 

ZXXl  =  CMPLX(X1,0.  0) 

ZXX2  =  (ALPA**2)*ZE  +  (BETA**2)*ZH 
ZXX  =  ZXX1*ZXX2 
X  =  ALPA*W/2. 

J31  =  2.*PI*SIN(X) 

J32  =  X**2  -  PI**2 

J3  =  J31/J32 

JXl  =  CMPLX(J3,0.  0) 

F4  =  ZXX*JX1**2 

RETURN 

END 
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APPENDIX  D.  PROGRAMS  FOR  CALCULATION  OF  FRINGING 


CAPACITANCE 


A.  TRANSMISSION  LINE  CIRCUIT  MODEL 

*  THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  THE  FRINGING  CAPACITANCE 

*  IN  SECOND  ORDER  APPROXIMATION  WITH  W  =  0. 7 11  MM  AND  ER  =  2.  2 

*  -  - -  -  ,  ■  - - - - -  ..  „  — 

COMPLEX  ZIN 

REAL  MAGZIN,OMG,ZO,BETA,AL,C1,C,FREQ,T,ZZ,PI,ZZIN 
PARAMETER  (  PI=3.  141592654,  Z0=94. 088,  AL=0.  005  ) 

PARAMETER  (  BETA  =  301.921744  ) 

*  NOTE  :  FILL  IN  PARAMETERS  BETA  AND  ZO  FROM  TABLE  1  OR  TABLE  2. 
DO  10  Cl  =  0.8058,0. 8062,1.  E-6 
C  =  C1* *1.E-14 

DO  15  FREQ  =  12.858241,12.858241 
OMG  =  2*PI*FREQ*1.E9 
T  =  TAN(BETA*AL) 

ZZ  =  OMG*C*ZO 

ZZIN  =  -l.*(l-ZZ*T)*ZO/(ZZ+T) 

ZIN  =  CMPLXC 0.0, ZZIN) 

MAGZIN  =  CABS((ZIN)) 

WRITE(6,100)  FREQ, C, MAGZIN 
100  FORMATC'  ' ,2X,F16.  12,2X,F25.  23,2X,F14.  12) 

15  CONTINUE 

10  CONTINUE 

END 


B.  END  EFFECT  EXTENSION  MODEL 

*  THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  THE  FRINGING  CAPACITANCE 

*  IN  SECOND  ORDER  APPROXIMATION  WITH  W  =  0. 711  MM  AND  ER  =  2.  2 


REAL  FO,BETA,LL,L,PI,DELTA,CO,OMG,ZO,CP 
PARAMETER  (  Z0=94. 088,  L=0.  01,  PI=3. 141592654  ) 

PARAMETER  (  F0=12.  858241 ,  BETA=301.  921744  ) 

*  NOTE  :  PARAMETERS  FO,  BETA,  AND  ZO  CAN  VARIES  DEPEND  ON  THE  CHOICE 

*  OF  APPROXIMATION 

OMG  =  2.*PI*F0*1.E9 
LL  =  PI /BETA 
DELTA  =  (LL-L)/2. 

CO  =  BETA/(OMG*ZO) 

CP  =  C0*DELTA 
WRITE(6  100)  FO, BETA, CP 
100  FORMATC'  ’ ,3X,F10. 7,5X,F9. 5,3X,F30. 29) 

END 
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APPENDIX  E.  PROGRAM  FOR  CALCULATION  OF 
CHARACTERISTIC  IMPEDANCE  OF  SUSPENDED  SUBSTRATE 

LINE 


THIS  IS  AN  EXAMPLE  PROGRAM  TO  COMPUTE  THE  DYNAMIC 
CHARACTERISTIC  IMPEDANCE  OF  SUSPENDED  SUBSTRATE  LINE 
IN  SECOND  ORDER  APPROXIMATION  WITH  W=l. 25  &  ER=9. 8 


REAL  EEFF,0MG,FREQ,BETA,W,G,ER,T,A,Z0,Z11,Z1,V,R 
PARAMETER  (  W=l. 25E-3,  G=1.575E-3,  ER=9.8,  T=0. 255E-3  ) 
PARAMETER  (  A=3. 2E-3  ,  BETA  =  293.61813  ,  PI=3.  141592654) 

*  NOTE  :  PARAMETER  A  IN  FIGURE  1  IS  2A  =  3. 2E-3  METER 

FREQ  =  9. 78036 

V  =  -1. 7866-0. 2035* *(T/G)+0.  4750*(A/G) 

R  =  1.  0835+0. 1007*(T/G)-0.  09457*(A/G) 

Zll  =  6*G/W  +  (1+4*((G/W)**2))**0.  5 
Z1  =  60*(V+R*L0G(Z11)) 

OMG  =  FREQ*PI*2. E+9 
EEFF  =  (BETA*3.  E+8/0MG)**2 
ZO  =  Z1/(EEFF**0.5) 

WRITE(6,25)  ZO 
25  F0RMAT(‘  ',5X,F9.3) 

STOP 

END 
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APPENDIX  F.  SOME  MATHEMATICAL  CALCULATION 


These  equations  are  used  in  computer  programs 
1.  The  equations  (Jp,q,m)  on  page  7 


>^1  =y«£o 

Zi  =j03H^ 

(1) 
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2.  The  equations  (ln,o)  on  page  7 
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3.  The  equation  (Ig)  on  page  6 

Where  : 

NE  =  yyiCt^  + 


7a  +  -  kI 

NE  =  - : - X  Ch  + 

1  Ct 


7?+7 


2  ^2 
^3 


X  Ch 


a£=^c-^7“'+/^'-^I  +  Cf27«' + /^' -  ^3  ] 
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X  y,  +  O,  X  Vs] 


(10) 


(11) 

(12) 

113) 

(14) 

(15) 

(16) 


4.  The  equation  (li)  on  page  6 

O]  =  Coth  y^h  =  Coth  ~  k]  x  h)  (17) 

IF  (a^  +  —  K^)  <  0,0  ,  then  y,  is  comple.x.  We  defined  y,  as  (7,,  so 
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(18) 


Oj  =  Coth  jG^h  = 
1 


1 


Cr,= 


Tank  jG^h 
-J 


JTan  G^h  Tan  G^h 


Oj  =  Comp/ex(0.0, 


ELSEIF 


C/j  =  Coth  G^h  = 


-1 


Tan  G^h 


1 


) 


lanh  G]h 


Cti  =  Complex{ 


1 


,0.0) 


Tank  <j|/z 

The  same  process  for  Ct^  and  Ct^  in  equation  (Ij)  and  (Ik). 

5.  The  equation  (Ih)  on  page  6 

NH 

Where  : 

NH  =  y22^h 


yl  +  /?^  —  K2  yj 

NH  =  ^ - X  C to  + 


+  f-  k] 


X  Ct-t 


jcop^  Mr  V 

We  assumed  that  the  value  of  p,  =  I.O,  then  we  have 
A  H  =  ~  ^Ctj  X  yo  Ct-i  X  ya  J 

jo)p^  ^  ^  ^  ^ 

Write  NH  in  a  symbolic  form, 

A’// =  ^-L_  X  [&&«&] 

J(^ho 


-  yzxyzi^h^h  ■+■  yz\yzi^h^h  +  yziyzz^h^h  +  v'r2 

/ 

V 


DH\  = 


/?+  -  k\  Ja^  +  P^-  kI 


yw/i. 


X  Ct^Ct2 


(19) 

(20) 

(21) 

(22) 


(23) 

(24) 

(25) 


(26) 

(27) 

(28) 

(29) 
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+  k\  +  f-  kI 

DH2  =  — - : - X  - - - X  OjOj 


Dm  = 


V 


yw/^o 

/ -  Kj  J<x^  +  (i^-Kj 


X  Ct2Ct^ 


DH4  =  [ 


Simplify  (  29  to  32  )  : 


D//1  = - ^  X  +  P^-Kj  X  x  0102 

{j^o) 


DH2  = - ^  X  +  /J^  -  kI  X  -  A'|  x  Ct^Ct^ 

Dm  = - X  +  /J^  -  X  +  /J^  -  A'J  X  Cr2Cr3 

O'^/^o) 

Z)//4  - - X  +  /J^  -  A'l  X  ^/ot^  +  /;^  -  a| 

(/WA^o) 


(31) 

(32) 

(33) 

(34) 

(35) 

(36) 


/)//  = - X  {DHl  +  DH2  +  Dm  +  DH4)  (37) 

iJ^CO^^o) 

Write  DH  in  a  symbolic  form, 

DH  = - X  tSSS]  (38) 

(/■w/'o) 

We  defined  in  the  program  that  Kf  =  =  AKl  and  Af  =  A K2 

^a^  +  P^-Kl  X  ^0^  +  P^-K]  =  AKK\ 

+  P^  -  AK^  =  AKK\  (39) 

and 

-  a|  X  -  A'|  =  /1AA2 

+  P^  -  AK2  =  AKK2  (40) 

From  these  equations  we  have  the  value  of  ZH  as  the  following; 
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z//  = 


X [&&&] 


[SSS] 


=yw^o 


[&&&] 

[sss] 


Defined  joJtx^  =  YYH 
YH  =  Complex{0.0,  YYH) 


So 


ZH=YHx 


[SSS] 


(41) 


(43) 
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